T • H • E 


OHIO 

SPOT 

UNIVERSITY 


STB-White, Final Report 


Student Authors: D. Molnar, R. Ammon, T. Gallagher, P. Gouhin, 
S. Hermann, J. Roos, C. Saurer, and H. White 

Supervisor: Dr. G. M. Gregorek 
Assistant: R. L. Reuss 

Department of Aeronautical & Astronautical Engineering 


Universities Space Research Association 
Houston, Texas 77058 

Subcontract Dated November 17, 1989 
Final Report 

RF Project 767919/722941 


May 1 990 



THE 


OHIO STATE UNIVERSITY 
AERONAUTICS DESIGN PROGRAM 


ADVANCED 


T • H • E 

OHIO 

SIAffi 


UNIVERSITY 



SPONSOR ED 


BY nasa/usra 



ilH 


STB — nri2± te 
Design Team 

Final Report 


Aero 694, 515H, 416 
May 30, 1990 



Dan Molnar 

- Group Leader 

— 


- Model Design 


Rob Ammon 

- Propulsion - Turbof anramjets 

__ 


- Landing Gear 


Todd Gallagher 

- Subsonic Aerodynamics 



- Supersonic Aerodynamics 


Pat Gouhin 

- Propulsion - SCRAMjets 

“ 


- Inlets 

v— 

Steve Hermann 

- Weight Analysis 

. _ 


- Stability and Control 



- Economic Analysis 


John Bryan Roos 

- Aircraft Design and Developm 

s 


- Waverider Design 

— 


- Environmental Aspects 


Craig Saurer 

- Materials 

— 


- Glide Analysis 

__ 

Heather White 

- Waverider Design 

r~ 


- Inlets 





ABSTRACT 


The final design of a hypersonic, SCRAMjet research aircraft, which 
is to be dropped from a carrier plane, is considered. Topics such 
as propulsion systems, aerodynamics, component weight analysis, and 
aircraft design with waverider analyses are stressed with smaller 
emphasis placed on aircraft systems such as cockpit design and 
landing gear configurations. Propulsion systems include analysis 
of the turbofanramjet for acceleration to low hypersonic speed 
(Mach 6.0) and analysis of the SCRAM jets themselves to carry the 
aircraft to Mach 10.0. Both analyses include the use of liquid 
hydrogen as fuel. Inlet design for both propulsion systems is 
analyzed as well. Aerodynamic properties are found using empirical 
and theoretical formulas for lift and drag on delta-wing aircraft. 
The aircraft design involves the integration of all preliminary 
studies into a modified waverider configuration. 
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1 . 0 T ntroduct ion 


For many years the concept of a supersonic combustion ramjet, or 
SCRAMjet, has been studied. SCRAMjets utilize high speed air flows 
to compress, burn with fuel, and expand air without the aid of 
moving parts such as compressors and turbines. The difficulty in 
designing and running these engines has prevented a SCRAMjet 
propelled aircraft from being built. Recent research has 

intensified, and the possibility of a working SCRAMjet is near. 
Work on the design of an aircraft to carry the SCRAMjet at the 
hypersonic speeds necessary for these engines to run is also 
intensifying. The purpose of this report is to present the final 
design work of a manned, SCRAMjet test aircraft, the STB-White 
(SCRAMjet Test Bed White Group). All aspects that have been 
analyzed for the last two quarters are included. Each section 
studied is included with its own discussion and conclusions. The 
main sections included are aircraft design, propulsion systems, 
aerodynamics and aircraft systems. The report also includes 
helpful figures such as 3-view drawings and pie charts at the end 
of each section to allow for simple study of data. 
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2.0 Aircraft Design and Developm ent 


2 . 1 Overview of Aircraft De s ign and Developmen t 

The STB-White was conceived as a SCRAMjet test bed. Its sole 

purpose is to carry a prototype SCRAMjet to a high speed, high 
altitude test run. This was kept in mind in all phases of design. 
All attempts were made to keep the aircraft small, light, and 
compact. The design of the aircraft completely centered around 
carrying the engines and the fuel to power them. Taking a look at 
Figure 2-1, it can be seen that the engines and the fuel take up 
60% of the total weight of the aircraft. The waverider 
configuration was seen solely as the most efficient way to get to 
the test run speed and altitude. The testing of a waverider 
configuration was seen only as a added bonus to the overall design. 

2 . 2 Mission Profile 

The mission requirements almost completely define the aircraft’s 
capabilities. The aircraft must be capable of: 

- dropped from a French designed carrier plane at Mach 0.8 

- climbing and accelerating under its own power to Mach 6.0 
and an altitude of approximately 80,000 feet 

- carry a prototype SCRAMjet engine that will start a Mach 6.0 
that will be able to climb and accelerate to Mach 10.0 and an 
altitude 100,000 feet, and maintain steady level flight for 
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two minutes at Mach 10. 

- decelerate and land at NASA Kennedy Space Center, Florida.* 

* This requirement was added by the STB— White design team to ease 
the problem of supersonic flight over land. The runway at the 
Kennedy Space Center is near the ocean. This will allow the flight 
to be completely subsonic over land. 

The aircraft meets these requirements. 

The current mission profile will be performed in five phases, 
transport, drop, acceleration, test run, and landing. Phase 1, 
transport, calls for the French carrier plane to transport the 
aircraft from Kennedy Space Center, Florida to Maine. In Maine, 
the aircraft will be fully fueled with LHg. Phase 2, drop, will 
have the French carrier plane flying to the drop zone approximately 
250 miles east and 100 miles north of Maine. Phase 3, 

acceleration, will include the acceleration up to Mach 6.0 using 
the turbof anramjets , and the acceleration to Mach 10.0 using the 
SCRAMjets. Looking at Figure 2-2, it can be seen that the SCRAMjet 
ignition is just south of Maine. In the case of the SCRAMjets not 
igniting the aircraft will be able to make an emergency landing in 
Maine. Phase 4, test run, will be steady level flight for two 
minutes at Mach 10.0 and an altitude of 100,000 feet. Phase 5, 
glide, has the aircraft slowing after the test run. It will 
approach the coast of Florida at around Mach 1.0, and continue an 
unpowered glide until it lands at Kennedy Space Center. 
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2 . 3 frpl imi nflrv Design 


The preliminary design, shown in Figure 2 3, was conceived 
the preliminary winter quarter presentation in mid- January. It was 
conceived based purely from concept. Few real calculations had yet 
been done, no SCRAMjet data was available at the time, little was 
known about waveriders, and the design was almost completely 
scrapped in favor of the next design. 


2 . 4 Hpsign Revisio n 1.0 


Design Revision 
the preliminary 
presentation, 
promise . 


1.0, shown in Figure 2-4, was a complete rework of 
design presented at the interim winter quarter 
This was the first design to show any merit or 


The main features of the aircraft is the over-under configuration 
of the propulsion system, and the waverider configuration. The 
turbofanramjets had been roughly sized, but little exact SCRAMjet 
data was available. Engine-out stability was the driving factor 
for the over-under configuration. Preliminary waverider analysis 
decided the basic shape of the aircraft. Optimization for Mach 6.0 
was decided upon because the high Mach number (8.5 - 10.0) data was 
not trusted, and the configuration remained relatively constant 

from Mach 6.0 - 8.0. 
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The position and type of the landing gear to be used was decided 
upon. The front gear would be located forward of the pilot and 
would be a steerable nose wheel. The rear gear, consisting 
pair of sleds, would be located behind the aircraft eg and outside 

of the SCRAMjets. 

A large problem arose from this design - internal space. There was 
not near the internal space needed to hold the LH 2 fuel. 


2 . 5 Design Revision 2.0 

Design Revision 2.0, shown in Figure 2-5, was presented at the 
final winter quarter presentation. This design solved the liquid 
hydrogen fuel volume problem. 


A complete rework of the waverider analysis had to be done. The 
length was extend 20 feet from the previous design to bring the 
total length to 80 feet. The span was also extended 22 feet for a 
total span of 56 feet. This allowed enough space for the LH 2 
needed for the mission. The second analysis of the waverider 
design produced more realistic data then the first, and it 
decided to optimize the aircraft for Mach 10.0. In this set of 
runs, the waverider configuration remained relatively constant 

after Mach 8.0. 


5 


The optimized waveriders from MAXWARP had somewhat pointed noses. 
It was felt that a pointed nose would be far to difficult to cool, 
and it was decided to stay with a blunt nose. Attempts were made 
to keep all changes in the planform smooth to aid in cooling. 

This design also shows the final configuration of the over-under 
propulsion system. It was decided that two downsized 
turbofanramjets and five full sized SCRAMjets could provide the 
necessary propulsion. 

2 . 6 Final Design Co nfiguration 

This design was formulated in the first few weeks of Spring 
Quarter. Over Spring break some flaws were found in Design 
Revision 2.0. It was almost impossible to fit the turbofanramjets 
and the SCRAMjets in the over-under configuration that was desired. 
Another problem was location of the cockpit. Variable geometry was 
needed for the turbofanramjets, and it would not allow any room for 
the cockpit. It was also decided to lengthen the aircraft to 85 
feet. This length was the maximum allowed by the French carrier 
plane. The current width was 56 feet, and the maximum allowed by 
the French carrier plane was 50 feet. The width was narrowed to 
48 feet for this reason. 

The waverider analysis was redone for the new length. The new 
geometry allowed enough room for the over-under engine 
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configuration, the fuel, and the cockpit. The final design was 
set, and final versions of all aspects of aircraft development 
began. The inlet study was begun, the model was designed and 
built, and stability and control were done. A 3-D model was 
created on I-DEAS V4 . 1 (Figure 2-6), which aided in the internal 
arrangement of component systems (Figure 2-7). 


Table 2-1 shows some of the more important specifications for the 
aircraft . 
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Specifications for the STB-White 
Table 2-1 


Length 85 feet 

Height 14 feet 

Span 4 ® feet 

Planform Area .2100 feet 

Wing Sweep ^0 degrees 

Aspect Ratio 1*1 

Wing loading 

At Drop... 25.2 1 bs / sqf t 

At Landing 18.6 lbs/sqft 

Gross DropWeight 53,0001bs 

Landing Weight 39,172 lbs 

Cruise Altitude 100,000 feet 

Cruise Mach Number 

Range 1400 miles 

Land ingVelocity 276 ft/sec (188 mph ) 

Landing Distance 5181 feet 

Propulsion 

Turbof anramjet 2 - 20% Downsized General Electric 

SCRAMjet 5 - Fullsized General Electric 
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2 . 7 WAVF.RIDER CONFIGURA TION 


When approaching hypersonic speeds, it is appropriate to have some 
kind of configuration which will yield as little drag as possible 
and, hence, more thrust available. For this reason, a partial 
waverider configuration was chosen for our design. Our analysis 
was done using the MAXWARP (Maryland AXisymmetric WAveRider 
Program) program which was setup on the Harris 800 at the Ohio 
State Aero/Astro Research Laboratory. 


In particular, the design was based on flow fields 
bodies. The following information was needed m an 
power-law bodies. 


over power-law 
input file for 


POWER 


In which, a power of the power law equation of 0.66 was 
used to obtain minimum drag characteristics. The power 
law equation is 


r 


■ base 


* (z/l)""' 


where rbase is the radius at z = 1, z is the axial 
distance from the nose of the power law body, 1 is the 
length of the body, power is the exponent of the power- 
law equation, and r is the local radtus of the power- 

1 aw body . 
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YBL = A value of 0.2 was the final number used for the base 

height/length ratio. 

ZPBODY = The length of the power law body was 18.29 meters (60 
f eet ) . 


ZSL = A value of 0.05 was used as the non-dimensional length 
of the nosetip to the length of the body ratio. 

Through '■'the iteration process, the length and base height were 
changed several times in order to assure that the volume was large 
enough to house fuel, engines, and aircraft systems. 


To account for different altitudes, the freestream conditions were 
changed to the appropriate values for various altitudes. The 
dynamic pressure and temperature were changed for several altitudes 
between 60,000 and 100,000 feet. Also, the specific heat ratio was 
changed between 1.36 and 1.4 to account for some disassociation at 
hypersonic speeds. It was found that changing this value had 
little affect on the design aspect of the configuration. Another 
factor in the freestream conditions was the Mach number. Changing 
this value did have a great affect on the design. For lower Mach 
numbers, the design was very rough, and the edges were very sharp. 
As the Mach number increased, the design became smoother, as shown 
in Figure 2-8. Also, the base height became thicker at lower 
speeds. This may be due to the increased size in the Mach cone on 
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the body. Our final configuration was taken from the Mach 10.0 
case from the waverider, and is shown in Figure 2-9. 

Aircraft analysis parameters were specified in the input file 
consisting of the aircraft length of 60 feet. Analysis was 
performed on a viscous body with turbulent air flow. The wall 
temperature was specified at 3000 degrees Rankine. 

The aircraft constraints were specified for volume and box size. 
A minimum volume was set at 5000 cubic feet, and the maximum was 
set at 7500 cubic feet. Box size is defined as the half span to 
length ratio. The minimum box size was set at 0.1, and the maximum 
box size was set at 0.4. These constraints were set to account for 
the amount of space that was needed for the aircraft fuel and 
systems. The slenderness ratio was set for a few runs, but the 
program kicked out these constraints. Therefore, for the final 
configuration, the slenderness ratio was not used. 

After coordinates were obtained from the H800 Super-mini Computer, 
these files were downloaded to a PC and read into an AutoCad V10.0. 
Once in AutoCad V10.0, the configuration was rotated and sized to 
obtain each view. 

MAXWARP was used to integrate several variables in order to obtain 
a design which minimized the drag count. The output conditions 
calculated from this program were only ideal values and could not 
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theoretically be used as final calculations. This was because the 
final configuration did not consist solely of the waverider 
configuration. Powerplants, boat tail, vertical stabilizer, and 
other details were added to the final configuration. So, the 
waverider was only used as a basis on which to design the 
hypersonic vehicle . 

2 . 8 Materials 


In order to chose materials, the thermal distribution is essential. 
The program CDHEAT (Conceptual Design aerodynamic HEATing analysis) 
finds this distribution. It uses algebraic equations based on 
empirical data and has a rather simplistic approach. However, for 
known body shapes and flight conditions it gives accurate results. 
Heat transfer coefficient ratios for the various body shapes are 
found and then used in an iterative process to find the heat 
transfer rates using the following equation: 


Q - h 0 -£- {T aw - T v ) 

h 0 - reference coef 
h - actual coef 
T mw - adibatic wall temp 
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Unfortunately, the code was not completed in time for this report 
but temperatures at the nose and along the leading edges were 
approximated. The nose temperature was under 3000 F and the 
leading edges were roughly 2400 F. This is beneficial since the 
aircraft probably will not need active cooling. This lack of an 
active cooling system will also save weight and space. 
Furthermore, since this vehicle experiences extreme temperatures 
for only a short time, the life of the surfaces at these 
temperatures need not be long. 

A carbon/carbon composite capable of withstanding temperatures 
over 3000 F for extended periods will sufficiently protect the 
nose. Its strength at high temperatures far surpasses any metal, 
and it can be used along the fuselage. At the leading edge, a 
refractory alloy of Columbium is effective. It is commercially 
available with reasonable strength at room temperature. It is also 
easy to fabricate and weld allowing easy attachment to our desired 
leading edge shape. Columbium oxidizes rapidly at temperatures 
over 1200 F. To prevent this, a coating of slurry silicide is 
fused onto the surface in a vacuum furnace. For the remaining 
"cool" spots, a titanium alloy can be used. 

Use of evacuated insulation for STB-Whi te ’s cryogenic fuel tanks 
was initially considered, but it was found that this type of tank 
often develops leaks in the insulation causing severe boil-off. 
The refore, a tank with internal insulation is superior. The system 
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consists of a heat shield surrounding a purge gap. The purge gap 
further surrounds the insulated tank on the inside. The only 
problem in this design is that LH 2 permeates the insulation raising 
the thermal conductivity and causing boil-off. To counter this 
effect, a vapor barrier is installed inside the insulation. 

2 . 9 Subsonic Wind Tunnel Model Construction 

At the beginning of Spring Quarter, Dr. Gregorek gave 
specifications for the subsonic wind tunnel model. A model two 
feet long made of wood (balsa or pine) that would mount on a 3/8 
inch stinger was requested. With these requirements in mind, a 
model was devised and constructed. 

The subsonic wind tunnel model was to be constructed of balsa wood. 

A 1/4" x 1/4" square spar was constructed. This spar would hold 
several rib cross-sections. Each rib was cut and placed in a 
direction perpendicular to the aircraft flight path - with the spar 
in the flight path direction. Next, 1/8" thick balsa sheets would 
be conformed to the ribs to make the skin of the aircraft. 

The front half of the aircraft was constructed first (i.e. the 
waverider portion). The underbody was simply formed by attaching 
the flat balsa sheets to the rib lower edges. The top of the 
aircraft was more difficult. The balsa sheet had to be bent to 
conform to the gently curving ribs that defined the aircraft upper 
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surface. To do this, water was added to the balsa sheets to make 
them more flexible. This flexible balsa wood was then wrapped over 
the rib contours and glued to complete the upper surface. 

The back half (i.e. the boat-tail) was completed in the same 
fashion as the front half. The only complication was adding the 
stinger cavity that must pass through the rear of the model for 
wind tunnel mounting. This cavity was created by cutting holes in 
the rib cross-sections and then lining them with sheet metal. 
Also, it is favorable to affix the stinger at the model center of 
gravity. This allows for simpler reduction of moment, lift, and 
drag data obtained from the wind tunnel tests. A hole was drilled 
through the model upper surface and through the sheet metal-lined 
cavity to allow a mounting point to the stinger. Next, the balsa 
sheets were conformed to the rib surfaces as before. Finally, the 
vertical control surfaces were attached to the model. 

Next, an attempt was made to make the model smooth. A micro foam 
filler was used to fill the balsa grain. The filler was sanded to 
smoothness with difficulties incurred at wood glue joints. The 
wood glue sanded more slowly than the balsa wood. The team logo 
was added to the model at this time by cutting it from a ’sticker’. 
The letters were glued to the model surface leaving the removable 
layer attached to the letters. The model was then coated with a 
clear coat enamel in an attempt to give the paint a less absorbent 
surface to adhere. Black enamel paint was added in several coats 
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and sanded and repainted to achieve a smooth surface. When 
painting was complete, the removable paper layer in the logo was 
removed leaving the logo embedded in the paint. 
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Fuel Tank 
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3 . 0 Propulsion Systems 


3 , i Overview of Propulsion Systems 


The propulsion system of the aircraft is broken down into two 
separate sections - turbof anram jets and SCRAMjets. The 
turbofanramjets operate from the drop at Mach 0.8 and an altitude 
of 40,000 feet to Mach 6.0 and an altitude of 80,000 feet. At this 
speed and altitude the SCRAMjets will ignite and continue 
accelerating and climbing up to Mach 10 at 100,000 feet. The 
aircraft will hold this speed and altitude for two minutes and then 
descend to earth. This flight profile is the most efficient 
integration of the two propulsion systems. 


3 . 2 Turbofanramjets 
3.2.1 Engine Analysis 

The engines that will be used to accelerate and climb the aircraft 
from Mach 0.8 at 40,000 feet to Mach 6.0 at 80,000 feet are Mach 
6.0 LHj-Fueled Augmented Turbofanramjet engines. These 
turbofanramjets will be produced by General Electric Aircraft 
Engines in Cincinnati, Ohio. The engines are scheduled to operate 
initially between the years 2005 and 2010. Dimensions, weights, 
and engine data for the Mach 6.0 engine were given to each of the 
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four design groups by General Electric. Two different Mach 6.0 
engines will be available, an acoustically unconstrained and an 
acoustically constrained one. A decision was made to use the 
acoustically unconstrained engine, because it weighed less and was 
smaller while producing the same thrust. Also, since our flight 
would take place away from land and people, the excess noise was 
not a factor in the design process. 

Last quarter, the flight profile (see Figure 3-1 and Figure 3-2) 
was determined and the thrust required and thrust available vs. 
Mach number plot (see Figure 3-3) was created. A decision was made 
to use two turbof anram jets for this portion of the mission. 
However, by using two full size engines, more thrust was available 
than was needed for this flight profile. Therefore, these engines 
could be sized down to cut down on weight and also give the 
aircraft more volume to store the liquid hydrogen fuel. General 
Electric included engine scaling relationships with the engine 
data. Upon examining the Thrust Required and Thrust Available vs. 
Mach number plot for two full size engines, it was determined that 
these engines could be down-sized by 20 percent with respect to 
thrust. The engine scaling relationships indicate that the thrust 
and fuel flow scaled one to one with the airflow. According to 
these formulas, the full size engine of 6100 lbs. scaled down to 
4880 lbs. (20%). Also, additional relationships for diameters and 
lengths were given and used. The results are as follows for a 20% 
reduction in engine size with respect to the thrust: 
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Table 3-1: Turbof anrarajet Size Reduction 


Full Size 20% Reduction 


Maximum Diameter (in.): 
Length ( in. ) : 

C.G. Location ( in. ) : 
Weight - each (lbs): 


83 

265 

119 

6100 


74 

237 

106 

4880 


The net thrusts given by O.E. can be reduced 20% 


Once the engines were scaled a configuration for the aircraft was 
made . The drag values that were finalised for our configuration 
were much higher than the preliminary ones of last quarter. These 
higher drag numbers gave the group two options. The first, an 
extra engine could be added to the aircraft to offset the enormous 
amount of drag in the transonic region. However, by adding an 
extra engine the entire aircraft would have to be sized up. This 
option was discarded because the French carrier aircraft design had 
given each group a size limit. If the aircraft was enlarged, the 
French limit would not be adhered to. Also, by changing the 
configuration of the aircraft, the drag analysis would have to be 
repeated and these drag coefficients would probably increase. 
Because of the uncertainty of the outcome in adding a third 
turbofanramjet engine, a second option was carried out. 


This second option was to alter the existing 
turbofanramjet portion of the mission, 
increased to almost double the available 


flight profile 


S ince 

the 

drag 

thrust 

of 

the 


for the 
values 
two 80% 
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turbofanramjet engines in the transonic range, the aircraft had to 
descend until the available thrust was appreciably greater than the 
drag. Once this was the case, the aircraft could ascend. Since 
the drag in the transonic region was so large, the aircraft had to 
descend to an altitude of 15,500 feet to overcome the drag and then 
ascend to 80,000 feet where the SCRAMjets were ignited. The 
updated flight profile and the Thrust Available and Thrust Required 
(drag) vs. Mach number are shown in Figures 3-4 and 3-5 and Figure 
3-6, respectively. 

3.2.2 Flight Profile 

The flight profile was analyzed by a more precise method than used 
in the past. The thrust available and drag values were known for 
the aircraft at different altitudes and at different mach numbers. 
The profile was altered by using only one parameter. This 
parameter, the angle of ascent or descent, was varied until 
positive acceleration was possible and when the fuel amount was as 
low as possible. The entire flight profile analysis was determined 
by summing the forces along and perpendicular to the flight path 
and setting these values equal to the mass times the acceleration 
(F=ma). The analysis was determined using a varying dynamic 
pressure (Q) until the aircraft ascended above 40,000 feet where 
the dynamic pressure was set to a constant of 1500 psf. A Q of 
1500 psf was determined to be the best for the profile. This Q was 
where the SCRAMjets performed well. Also this constant Q analysis 
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determined the velocity and thus Mach number at each station of the 
ciimb. The fuel estimate for this profile was calculated by 
knowing the thrust available (used), the specific fuel consumption 
Ufc) used by the engine, and the duration at which the engine used 

this fuel flow rate. Even though the previous engine data could be 

■ ei j A un i t was unfortunate that 
reduced when the engines were sized down, it was un 

the sfc could not also be reduced. If it could have been, the 

lowered sfc would have translated into a lower amount of fuel 

, • Thf* fuel needed to accelerate from Mach 

needed for the mission. The fuel neeaeu 

o 8 at 40,000 feet to Mach 6.0 at 80,000 feet was 6460 lbs of 

li q uid hydrogen. This amount of fuel translates into a storage 

volume of 1452 cubic feet at a density of 4.43 lbs per cubic foot, 
specifications of the flight profile are included in Table 3-2. 


Table 3 - 2 : Turbof anramjet Acceleration Data 


MAXIMUM ACCELERATION: 

TIME (PORTION OF TOTAL MISSION). 
TOTAL FUEL WEIGHT: 

AIRCRAFT T.O. WEIGHT: 

DISTANCE TRAVELED (PORTION): 


1.35 G 
6.10 MINUTES 
6,460 LBS 
53,000 LBS 
158 MILES 
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3 . 3 SCRAM.iets 


3.3.1 Approach To Analysis 

Above Mach 6, the only hope for an air-breathing engine that exists 
today is the operation of the SCRAMjet engine. At the present 
time, this idea is still in its preliminary stages.. Due to this 
fact, the data that is available on these engines is sparse and 
hard to acquire. The mission requirements of this design project 
call for the acceleration of a hypersonic vehicle to ten times the 
speed of sound at an altitude of 100,000 feet. Once at this speed 
and altitude, the vehicle must maintain steady level flight for two 
minutes. The airplane is to weigh approximately 50,000 lbs. and 
have an internal volume of 3,000 cubic feet allowed for fuel. With 
this brief summary of objectives, the possibilities will be 
analyzed, and the steps taken will be discussed. 

General Electric kindly provided the advanced design class with 
SCRAMjet data that was developed from recent tests. Without this 
data the following efforts could not have been performed and the 
whole mission would have been impossible. This data consisted of 
net thrusts that were developed by a liquid hydrogen burning 
engine. Since these thrusts were net thrusts, the ram drag was 
already accounted for and is not considered in this analysis. G.E. 
provided this information for dynamic pressures ranging from 500 
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psf to 2,000 psf in 500 psf increments. This covered altitudes 
from 65,000 feet to 126,000 feet. Combustor exit fuel to air 
equivalence ratios of 0.8, 1.8, and 2.8 were used in determining 
the specific impulse for different thrusts. The data covered Mach 
numbers of 5.0, 7.5, and 10.0. This fit into the regime of flight 
very well. Combustor exit static pressures and static temperatures 
were also given. The only regret was that more points in this 
range were not available for analysis. 

This data, along with information from different members of the 
group, allowed for calculations to be made. The aerodynamicist 
developed initial values for drag coefficients that were used to 
calculate the acceleration of the aircraft. Weights were estimated 
so the required thrusts could be calculated. From the waverider 
analysis, initial estimates were made on lift to drag ratios at 
different Mach numbers. The design procedure was a cyclical 
process that involved many iterations and refinements. 

3.3.2 Strategic Requirement Analysis 

From the G.E. data, the most thrust was developed for a dynamic 
pressure of 2,000 psf. This was the data that was used for the 
first round of calculations. However, it was soon realized that at 
this dynamic pressure with a velocity of Mach 10.0, the altitude 
was not as high as what was needed. From simple calculations, it 
was found that the dynamic pressure that would correspond to Mach 
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10.0 at 100,000 feet on a standard atmospheric day was 1,615 psf. 
The closest available data that could be used was at 1,500 psf. 
Since the turbofanramjet engines are operable to Mach 6.0, this 
point was decided on to ignite the supersonic combustion ramjet 
engines. Back solving the dynamic pressure equation showed at Mach 

6.0 the aircraft should be flying at approximately 80,000 feet. 
This meant the vehicle would accelerate from Mach 6.0 to Mach 10.0 
while climbing from 80,000 feet to 100,000 feet. The data from 
G.E. was modified so an estimate could be made for values of Mach 

6.0 at 80,000 feet and Mach 9.0 at 97,000 feet. Information was 
now available at Mach numbers of 6.0, 7.5, 9.0, and 10.0. 

3.3.3 Iteration Procedures 

Only four points were available in this region, so four divisions 
were made for calculations. The information from Mach 6.0 would be 
used from Mach 6.0-6.75. The information from Mach 7.5 would be 
used from 6.75-8.25. Mach 8.25-9.25 would use data from Mach 9.0, 
and finally Mach 9.25-10.0 would use the Mach 10.0 data. The 
initial weight of the airplane, at SCRAMjet ignition, was estimated 
at 50,000 lbs. All of the data was put into a spreadsheet and 
examined for the cases of using 1-8 modules (engines). It was 
quickly narrowed down to examine the realistic choices of 3 to 6 
engines. To see this data and the process that was involved refer 
to Appendix A. Data was analyzed concerning the following items: 
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Acceleration 

( # Eng x Net thrust - C D x Q x S ) / mass 
Time of operation 


( V final _ V initial 


) / Acceleration 


Fuel Weight 

( # Eng x Net thrust / I gp x 3600 x Time(in hours) ) 


By examining these three major factors, it was seen that the 
smaller the thrust, the smaller the acceleration. The smaller the 
acceleration, the longer the operational time. The longer the 
operational time, the more fuel that was used. There are 
exceptions depending on equivalence ratio and specific impulse, but 
this was the general trend. The two minute steady level flight was 
not considered in these calculations. This was because at cruise, 
the objective was only to offset the drag. This was possible to do 
with only two engines. However, by shutting down three of the 
SCRAMjets, the drag would greatly increase as no flow was directed 
through the engines. It was quickly realized that running all five 
engines at a low fuel to air mixture equivalence ratio would solve 
the problem. No data was supplied from G.E. for equivalence 
ratios below 0.8. Approximate data was generated for lower values 
to get a rough estimation that 600 lbs of fuel would be burned to 
offset the drag if all five engines were operating. 
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3*3.4 Discussion of Results 


It was realized that operating with less than five engines would 
consume an enormous amount of fuel and operating above five engines 
would not save enough fuel to offset the 1200 pounds that another 
SCRAMjet module would weigh. The next problem that was faced was 
an overabundance of lift. This made it impossible to fly at the 
requested flight profile. In order to overcome this problem, the 
angle of attack was changed until an achievable flight profile was 
reached. This occurred at an angle of attack of -4 degrees. 
Running through the spreadsheet with the corresponding coefficients 
of drag for this angle of attack yielded the final results. It was 
decided to use an extra 500 lbs. of fuel over the optimized value. 
This was done to increase acceleration and shorten range and flight 
time. By upping the fuel to air mixture equivalence ratio from 0.8 
to 1.8 when the SCRAMjets are accelerating from Mach 6.0 to Mach 
8.25, the range was shortened by 175 miles and the total flight 
time was decreased by 2.13 minutes. To accelerate from Mach 8.25 
to Mach 10.0 the optimized value was still used. which occurred at 
a ratio of 1.8. The total range over which the SCRAMjets operated 
was 584 miles. The vehicle took 5.82 minutes to travel this 
distance and used 7428 lbs. of LHj in the process. All of these 
numbers take into account the two minute cruise. 



Mach Number vs. Time (Last Quarter) 


gure 3- I 



joqurnH 
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Thrust Req. & Avail, vs. Mach No. 
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Mach Number 




Figure 3-4 
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Thrust Required, Thrust Avail, vs Mach 
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4 . 0 Inlet Design 


4 . 1 Overview of Inlet De sign 

Inlets are a series of ramps which direct the air flow into an 
engine. Axisymmetric inlets usually consist of flow over a cone in 
which ramps are formed from the basic cone shape, whereas, two- 
dimensional inlets consist of flow over wedges. Generally, two- 
dimensional inlets are considerably easier to analyze. The purpose 
of an inlet is to deliver the flow to the engine such that the flow 
is perpendicular to the engine compressor face. Therefore, the 
ramps are devised to turn the air flow to meet this condition. 
However, this is not the only function of the ramp system. The 
shocks that are caused by the ramps are needed to slow the flow 
down before entering the engine. Unfortunately, mass flow rates 
cannot be matched at all times during flight without having an 
extravagant inlet system which consists of many moving parts. The 
number of moving parts is directly related to the structural 
weight, so it is impossible to have this "flexible wave" which can 
be transformed into any shape. On the other hand, if the inlet 
capture area is not matched up with the needs of the engine, many 
problems could occur. If not enough flow is delivered to the 
engine, then the engine will not operate efficiently; if too much 
flow is delivered to the engine, the air flow will build up and 
spill out causing spillage drag. There are some ways to correct 
for the two air flow problems by adding air ducts and vents. 


42 


Besides all the described problems, boundary layer build up must 
also be considered when designing an inlet. Generally, the 
boundary layer is bled off with a splitter plate which is placed a 
prescribed length parallel to the body. This allows a new growth 
of boundary layer for a smaller distance on the splitter plate. 
The losses due to the newly formed boundary layer are not as great 
as ~ the boundary layer formed from the length of the aircraft. 
Also, boundary layer suction can be performed using a porous plate 
in which the change in pressure will "suck" the boundary layer off 
the surface . 


The main objective of the test vehicle was to test a SCRAMjet 
engine for two minutes at Mach 10.0 at an altitude of 100,000 feet. 
After being dropped from a French designed carrier aircraft, a 
turbofanramjet engine will be used to propel the aircraft from 
40,000 feet at Mach 0.8 to 80,000 feet at Mach 6.0. Since the two 
80 percent full size turbofanramjet engines are placed inside the 
vehicle, the flow needs to be turned into the engine, whereas, the 
five full size SCRAMjet engines are located on the lower surface of 
the STB-White. From the design of the plane, all the engines are 
sharing one inlet surface on the lower side of the plane. Using a 
series of doors and ramps, the complicated system will have the 
turbofanramjet engine closed off during the SCRAMjet test. Since 
the STB will be flying at a negative angle of attack of -4 degrees 
during the entire time the SCRAMjets will be operating, they will 
be angled 4 degrees to the top surface of the plane. This will 
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For 


allow for higher pressure recovery for the SCRAMjet inlet 
both the turbofanramjet and the SCRAMjet inlet, a complicated 
program was written in order to obtain conditions at the entrance 
of each engine. Static pressure and temperature, combustion Mach 
number, total temperature and pressure, capture areas, and mass 
flow rates were calculated from freestream density, pressure, 
temperature, and Mach number which were input along with the 
desired ramp angles. Following this analysis of each inlet, an 
appropriate result was obtained for both engine inlets. (see 
Appendix B for turbofanramjet inlet Fortran code and Appendix C for 
SCRAMjet inlet Fortran code). 

Because of the basic wedge shape of the waverider design, a two- 
dimensional inlet using linear compression ramps and variable 
geometry will be used for both the turbofanramjet and the SCRAMjet 
inlets. Linear ramps are used for their ease in analysis. 


4 . 2 Turbofanramjet Inlet 


Three constraints were placed on the turbofanramjet inlet design 
based on conditions at the compressor face. First, the airflow 
needed to be slowed down to subsonic speeds before entering the 
compressor. Second, the capture area at the compressor face had to 
tie a 3 feet radius for each engine, and third, the flow conditions 
had to match the engine data provided by General Electric for 

various speeds and altitudes. 
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Two design Mach numbers were chosen for the design of the 
turbof anramjet . First, Mach 2.0 was analyzed using the computer 
code. Then Mach 4.0 was analyzed in the same manner. The program 
used oblique shock relations and normal shock relations obtained 
from Ref. 2 to analyze a series of ramps and angles. First, the 



satisfactory result was obtained. The number of ramps was set at 
three due to geometry constraints on the lower surface. 


With the flow conditions known for the engine, guesses were made 
for ramp angles for both design Mach numbers. For the Mach 2.0 
design condition, by using an initial cone angle of 3.4 degrees (8 
degree wedge flying at -4.6 degree angle of attack), an initial 
ramp angle of 5 degrees, and cowl angles of 7 degrees each, the 
speed was slowed down to Mach 0.838 with a pressure recovery of 
0.981. For the Mach 4 condition, the initial cone angle was set to 
5 degrees, the initial ramp angle was 14 degrees, and the cowl 
angles were 16 and 18 degrees. This slowed the air flow down to 
Mach 0.766 with a total pressure recovery of 0.972 (see Appendix 
B). As one can tell, the construction of the lower surface of the 
STB is going to be quite complicated with variable geometry 
compression ramps. Following the calculation of the compression 
angles, the distances from the nose of the aircraft were calculated 
where the ramps needed to be placed in order to impinge on the cowl 
lip to form the compression inlet. In order to turn the flow up 
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into the turbofanramjet, the flow needed to be subsonic prior to 
the turn. Therefore, a normal shock formed on the cowl lip (see 
Figure 4-1). This caused an external compression that was not 
desired, but was considered to be less of a problem than the 
expansion that would have otherwise resulted. This expansion would 
have caused the flow velocity to increase, and the normal shock at 
the compressor face would have caused entropy gradients that were 
extremely high. Since the mass flow rate did not match up for the 
engine, a number of vents were placed on the engine to route the 
excess air around the entrance of the compressor. 

In conclusion, the turbofanramjet inlet did not exactly match the 
engine needs as far as mass flow rate conditions. It must be noted 
that the General Electric engine data was unclear as to the 
definition of. the engine mass flow rate. As to which is more 
costly to the design regarding weight penalties and spillage drag, 
the answer is beyond the scope of this report. Not enough 
knowledge is known pertaining to the weight of variable geometry 
inlets. Therefore, this design entails many losses due to weight 

penalties . 

4 . 3 SCRAM i e t Inlet 

Like the turbofanramjet inlet, the SCRAMjet inlet has certain 
constraints which need to be met. First, the air flow to the 
engine needs to be slowed down to between Mach 3 and Mach 5. 
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Second, the capture areas need to be matched to the full-size 
SCRAMjet, and third, the mass flow rate of the air needs to be 
equal to that of the required air mass flow rate of the engine. 

Again, a linear ramp system with variable geometry is used for the 
SCRAMjet. With the aid of the computer code, the pressure 
recovery, mass flow rate of the air flow, and the speed of the air 
flow to the engine were calculated using freestream conditions and 
ramp angles. The oblique shock angles were calculated using 
supersonic oblique shock relations. These were the simplest and 
best approximations available to calculate the shock angles. The 
angles were iterated using the theta-beta-Mach number relation (see 
Ref. 2). 

The SCRAMjet inlet was designed to handle many speeds in order to 
insure fewer losses due to spillage drag. Since the primary 
function of this aircraft was to test the SCRAMjets, the inlets 
were designed to eliminate as many air flow problems as possible. 
The design Mach numbers were chosen as Mach 6.0, Mach 7.0, Mach 
8.0, Mach 9.0, and Mach 10.0. A three ramp system was chosen for 
the SCRAMjet. The turbof anram jet inlet will be completely closed 
off by connecting the cowl and the inlet surface. This will give 
a virtually smooth surface to work with for the SCRAMjet. The 
shock which occurs off the nose of the STB-k*hite will completely 
miss the bottom edge of the SCRAMjet at all times. Therefore, the 
ramp system begins with one ramp on the lower surface of the 
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aircraf t . 


The other two ramps are located on 
SCRAM jet. Three ramps were chosen on account 

geometry system limitations. 


the cowl of the 
of the variable 


For the Mach 6.0 design case, the cone angle was set at 4 degrees 
(8 degree wedge flying at -4 degree angle of attack), and the first 
ramp angle was 6.5 degrees. The two cowl angles were 7 and 3.5 
degrees. This gave a total pressure recovery of 0.885 and a speed 
of Mach 3.8 entering the SCRAM jet. The pressure inside the 

combustion chamber was 703.0 psf , and the mass flow rate was 
matched at 142.0 pounds per second. Therefore, the intake capture 
area was matched to the SCRAMjet inlet area at 8 inches high and 30 
inches wide. For the Mach 7.0 case, the cone angle was 4 degrees 
with the first ramp angle at 5 degrees. The cowl angles were 5 and 
4 degrees respectively. This resulted in a total pressure recovery 
of 0.90 and a speed of Mach 4.5. The pressure inside the 

combustion chamber was 560 psf, and the mass flow rate was matched 
at 135 lb/s. For the Mach 8.0 case, the first cone angle was 4 

degrees, and the first ramp angle was 5 degrees. The cowl angles 

were calculated to be 6 and 3 degrees respectively. For the Mach 
9.0 case, the first cone angle was 4 degrees with a first ramp 
angle of 4 degrees. The cowl angles were both 4 degrees. For the 
Mach 10.0 case, the first cone angle was 4 degrees with a ramp 
angle of 3.5 degrees. The cowl angles were 3.5 and 2 degrees, 
respectively. See Appendix D for the computer code output that is 
referred to above. 
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In conclusion, the turbof anram jet engine took much larger losses xn 
spillage drag than the SCRAM jet. The SCRAM jet flow was matched for 
every design Mach number. This was done in order to concentrate on 
the abilities of the SCRAMjet without complications in air flow. 
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5 . 0 Aerodynamics 


5 . 1 Subsonic /Super sonic /Hyper sonic Aerodynamics 






f=I 


Working independently from the engine groups and hoping to cover 
any future changes in flight profile the approach to the 
aerodynamic problem was formulated for a range of possible Mach 
numbers and angles of attack at varying altitudes. Dimensionless 
aerodynamic parameters were found through shock expansion theory , 
DATCOM methods, and Nicolai’s approximations. Comparison of the 
different methods showed that the DATCOM 2-D wedge approximations 
resulted in the highest C La ’ s , The comparison also showed that 
Nicolai’s and DATCOM ’ s delta-wing approximations corresponded to 
one another almost exactly. The delta-wing C La’ s were sign if icantly 
lower than those from the 2-D wedge, however, it was believed that 
they represented the actual flow field the best. Therefore, the 
delta-wing approximation was used to calculate the supersonic and 
hypersonic aerodynamic characteristics of the STB-White . Subsonic 
characteristics were approximated assuming a linear lift curve 
slope . 


Our goal was accomplished by using a spreadsheet, which calculated 
the C L ’ s , Cj’s, and L/D’s for every Mach number at each altitude and 
angle of attack. Subsonic and Supersonic C L ’s were calculated from 
equations B.l and B.2 respectively. Subsonic ’ s were calculated 
using the aspect ratio, Mach number, and sweep angle. The 
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supersonic and hyperson ic C La’ s came from a delta-wing lift curve 
graph provided by Nicolai and DATCOM. -The C^ & ’ s were then 
incorporated with angles of attack ranging from —4 to 6 degrees, 
to calculate a variety of ’ s at each Mach number and altitude. 
From the density and Mach number at each altitude, a velocity and 
Reynolds number were calculated. Using equations B.3 and B.4 
subsonic, and super /hypersonic C Dq 1 s were calculated respectively. 
Subsonic Cp Q 1 s were based exclusively on friction drag, having 
Reynolds number as their only parameter. Supersonic and 
hypersonic Cp Q 1 s were composed of friction drag and wave drag , 
having Mach number, thickness to chord ratio, Reynolds number, and 
the B factor as their parameters. With C Do 9 s and C L ’s for each Mach 
number, angle of attack, and altitude, the total Cp was calculated 
using equation B.5. The total Cp’s were then multiplied by an 
efficiency factor of 1.1 to account for cowl drag along with any 
other miscellaneous drag items. The aerodynamic calculations done 
here considered the entire upper surface of the STB-White to be a 
lifting surface. The expansion lift from the exhaust ramp was 
neglected. Inquiries were made into the effects of centripetal 
lift, however, they were small enough to be considered negligible. 
Aerodynamic graphs can be seen on the accompanying pages. To see 
data, refer to Appendix F. 
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5 . 2 Glide Analysis 


An unpowered glide and land was designated after the two minute 
SCRAMjet test to cut down on fuel volume. The glide trajectory was 
calculated with a spreadsheet using aerodynamic data. 

The spread sheet contained the altitudes from 100,000 feet to 0 
feet in 1000 feet increments and their corresponding densities. 
The analysis began with a dynamic pressure, Q, of 1500 psf at 
100,000 feet and was gradually lowered to 60 psf at 0 ft. From the 
dynamic pressures and the density, a velocity could be found. A 
free body diagram gave the equations to solve for lift and 
acceleration. The angle of descent was arbitrarily chosen. 


L = W cos 0 


ma = W sin 0 -D 

Cl = W cos 0 / q S 

Next, the dimensionless lift and drag were obtained from our 
ae rodynam i c i s t . The CL and CD values gave the angle of attack 
using the equation: CL = CLo (a - ). Then, acceleration, or 

in our case, deceleration, was found using the dynamic equation 
F = ma . Now, using the equations of motion: 
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' ff" 


vf ~2 - Vi"2 .= 2as 


s 


= Vi t + .5 a t~ 2 


where s - distance travelled 


Vf = final velocity 


Vi = initial velocity 


The distance travelled and the expended time were calculated. This 
distance was along the path travelled, so the x-y components could 
be found using the angle of descent. 


Sx = s cos e 


Sy = s sin © 

Now, this distance in the y-direction should equal 1000 feet. When 
it did not equal 1000 feet, the angle of descent was varied until 
: it did. These calculations were made for all cases. 
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CL vs MACH 



MACH NUMBER 







□ c 





CD vs MACH 

Alpha*** —2, O, and 2 at 60,000 ft 
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MACH NUMBER 
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000 ft 



L/D vs ALPHA 

Mach= .6 at 40K, 60K, 80K, and 100K ft 


I./D vs, A ]]:ha - Mach 0.6 
Figure 5-5 
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L/D vs ALPHA 

Mach® 4 at 40K, 60K, 80K, and 100K ft 


L/D vs. Alpha - Mach 1.0 
Figure 5-6 
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. ALPHA (degrees) 

40,000 ft -+- 60,000 ft 80,000 ft -B- 100.000 ft 



L/D vs ALPHA 

Mach= 6 at 40K, 60K, 80K, and 100K ft 


- Mach . 0 
Figure 5-7 
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ALPHA (degrees) 

40,000 ft -+-60,000 ft 80,000 ft 100,000 ft 






CD vs CL 

Mach= .6 at 40K, 60K, 80K, and 100K ft 



40,000 ft — I — 60,000 f t -*“80,000 ft -a- 100,000 ft 





CD vs CL 

Mach= 2 at 40K, 60K, 80K, and 100K ft 



40,000 ft -*-60.000 ft -*-80,000 ft' 100.000 ft 
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6 . 0 Stability and Contro l Analysis 


6 . 1 Weight Analysis 

A preliminary estimate of the take off weight of the aircraft was 
made using historical trends for past experimental and high speed 
aircraft. This method was outlined in Reference 9. This 
preliminary value for the take-off weight was then used in 
subsequent configuration calculations. After a configuration was 
determined, it was possible to determine a more accurate take-off 
weight for the aircraft using a component breakdown. A new 
configuration was then calculated using the value for the take off 
weight determined by the component breakdown method. This process 
was repeated several times to produce an accurate result. The 
final value of the take off weight was 53,000 lbs. 


The basic design of the aircraft was determined by the waverider 
program plots. Given the basic outline of the aircraft, it was 
then necessary to configure the aircraft for optimum stability and 
fuel volume. A preliminary design proved to be inadequate for 
meeting the fuel requirement, so the entire design was enlarged. 
The final configuration consisted of a front mounted nose wheel, 
followed by the cockpit, with the twin t u rbo f anram j e t s placed side 
by side in the rear. Conformal fuel tanks were placed in the wing 
area of the aircraft. In addition to meeting our fuel volume 
requirements, the final configuration also included several mission 
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considerations . 


The turbofanramjets were mounted side by side to 
increase engine out stability and to allow one inlet to be used for 
both turbofanramjet engines. An inlet system length of eleven feet 
was accounted for in the configuration. The cockpit dimensions 
included adequate space for the data gathering computers and the 
cockpit cooling systems. In addition, a coverable window was 
included for use during landing. The nose wheel was mounted with 
a forward hinge to permit the wheel to be in proximity with the 
cockpit cooling system. The placement of the fuel tanks allowed 
for a small center of gravity (c.g.) travel. This was because the 
fuel tanks had a c.g. which was very close to the c.g. of the 
structure . 

Once the final configuration was determined, it was possible to 
perform a component breakdown to determine the weight of the 
aircraft. The basic areas of weight on the aircraft were; 
structure, engines, fuel, and systems. The fuel weight was 
determined by the amount of fuel required for the mission ( 13960 
lbs), and the basic weights for the engines were determined to be 
10000 lbs for the ramjets, and 6000 lbs for the SCRAMjets. 
Estimates for the weight of the remainder of the components were 
made using the weight estimation methods outlined Reference 9. 
These estimates were based on parameters involving the geometry, 
aerodynamics and estimated take-off weight of the aircraft. The 
weight values for the inlet, vertical tails, computer systems and 
cooling systems were estimates made from historical trends, 
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combined with information determined in our studies of these 
components. Pie graphs of the weight breakdowns for total weight, 
structural weight and systems weight are shown in Figures 6-1 
through 6-3 . 

Once the component weights were determined, it was possible to 
determine the c.g. of the empty plane by performing a moment 
balance on the plane. This involved adding the product of the mass 
and c.g. distance of each component, and then dividing this sum by 
the total mass of the aircraft. The empty c.g. of the aircraft was 
determined to be 53.2 feet from the nose of the aircraft. The c.g. 
of the full plane was determined by also adding in the fuel 
mass-c.g. product into the c.g. determination. The c.g. of the 
full aircraft was found to be 51.4 feet from the nose of the 
aircraft. Thus, the total distance traveled by the c.g. was only 
1.8 f ee t . 


The component breakdown of the weight of the aircraft, and the c.g. 
determination of the aircraft are presented in Appendix H. 
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6 . 2 Stability and Control 


6.2.1 Overview of Stability and Control 

In determining the stability of the aircraft, the main concern was 
for static stability at low speeds, particularly landing speeds. 
The main reason for this concern was due to the fact that the 
aerodynamic center was at it’s most forward position at low speeds, 
which produced the smallest static margin, and therefore the least 
stability. Also, the stability of the aircraft would be critical 
during landing to ensure a safe return for the aircraft and the 
pilot. Dynamic stability was not checked, but stability 

augmentation systems will be used to correct dynamic instabilities. 
The static stability of the aircraft was determined using the 
methods presented in Reference 9. The various stability 

derivatives were determined using References 9, 10, and the USAF 

DATCOM volumes on aircraft stability. Since the aircraft had no 
horizontal tail, elevons were chosen to provide both longitudinal 
and lateral control. For directional control a pair of vertical 
stabilizers were mounted in the rear of the aircraft. For this 
configuration, it was necessary to use the stability equations for 
a tailless aircraft. The aircraft was checked for stability in all 
three planes of motion: longitudinal, lateral and directional. The 
stability of each of these planes will be discussed separately. 
The aircraft was found to be stable in all three planes, and some 
performance criteria were determined. A listing of the important 
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static stability parameters is shown in Table 6-1, and a complete 
listing of the stability and control calculations is provided in 
Appendix I . 


6.2.2 Longitudinal Stability 

The primary factor in determining longitudinal static stability was 
the distance between the aerodynamic center and the center of 
gravity, or the static margin. For a tailless aircraft, the 
aerodynamic center and the center of gravity are the same point, 
and for longitudinal stability it is essential that the center of 
gravity, or c.g. , be in front of the aerodynamic center. However, 
if the static margin were to become too large, the aircraft would 
be too stable, and thus hard to control. To prevent this it was 
necessary to ensure that the static margin remain relatively 
constant throughout the flight. This will be accomplished by fuel 
staging to match the flight profile of the aircraft. 

The positive longitudinal stability criterion is that the slope of 
C ^ vs. angle of attack, or , be negative. This insures that the 
aircraft will tend to correct its pitch attitude when disturbed 
from the equilibrium position. The value of determined for the 
aircraft was -0.105, thus meeting this condition. Also, the value 
of Cttrr at zero angle of attack should be positive to ensure 
longitudinal stability. The value of at zero angle of attack 
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was determined to be 0.03> and this condition was met also. A plot 
°f C„ CG vs. angle of attack is shown in Figure 6-4. 

6.2.3 Lateral Stability 

The primary requirement for lateral stability was that the lateral 
stability derivative, or Cy , be negative. This requirement not 
only insures positive lateral static stability, but also keeps the 
spiral mode from becoming divergent. The value of Cy determined 
for the aircraft was —0.726, thus meeting this condition. For 
calculation of the steady state roll the primary criterion was the 
damping in roll coefficient, or C Lp . A value of -0.105 was 

determined for Cy , thus producing a roll rate of 87 deg/sec upon 

landing. This high roll rate was due to the low aspect ratio of 

the aircraft, and also because the elevons were sized for 

longitudinal stability, which provided a great deal more lateral 
control that would otherwise be present. A plot of roll rate vs 
velocity is shown in Figure 6-5. 

6.2.4 Directional Stability 

The primary requirement for directional (weathercock) stability was 
that the directional stability derivative, or , be positive. 
This requirement insures that the moments generated will rotate the 
aircraft to reduce the sideslip angle, or D. Cj| t was determined to 
be 0.102, thus this condition was met, also. 
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For sizing the aircraft rudders, it was necessary to check for 
adequate control under two conditions: maintain B of zero for one 
engine out, and maintain 13 of zero for landing in a crosswind. 
From this, a maximum value of x was determined, and the rudder area 
could be calculated. The total rudder area was about 125 square 
feet, or 0.28 of the total vertical stabilizer area. 

6.2.5 Trim Requirements 

The trim requirements were determined from the elevon control 
power, or ECP. The required elevator deflections for flight at 
various angles of attack could be determined by use of the ECP. A 
plot of the required elevator deflection vs angle of attack is 
shown in Figure 6-6. 

The turn radius required for various maneuvers was also calculated 
using force balance equations in combination with the control 
parameters. The maneuvers were the turn, pull-up and pull-down. 
The maneuvers were calculated for an artificial gravity of 3 g s. 
Plots of the radii vs velocity are shown in Figure 6-7. 
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Table 6-1 


STATIC STABILITY PARAMETERS 

LONGITUDINAL STABILITY 
Static Margin (Fully Fueled): 
Static Margin (Empty): 

CMA (Slope of CM vs Alpha): 

CM at Zero Angle of Attack: 

LATERAL STABILITY 

Rolling Moment Coefficient: 

Lateral Stability Derivative: 


NEEDS TO BE +/- 
FOR STABILITY 


0.139 + 

0.094 + 

-0.105 

0.030 + 

0.069 + 


- 0.726 


DIRECTIONAL STABILITY 
Directional Stability Derivative: 


6 . 3 Subsonic Wind Tunnel Testing 

Subsonic wind tunnel testing was to be completed using the subsonic 
wind tunnel model. When mounting was attempted, certain problems 
were found. First, the model was too long and the stinger was not 
long enough to reach the model center of gravity. Next, the model 
has a large flat plate lifting surface. This would be a problem 
when the model was placed at high angle attack possibly destroying 
the delicate strain gages under high loads. With these problems, 
the model needed to be torn apart and a different mounting 
technique used. Since there was little time left in the project, 
these modifications could not be made and the subsonic wind tunnel 
testing could not be completed. 
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Roll Rate vs. Velocity 
Figure 6-5 
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7 . 0 Aircraft Systems 


7 . 1 Cockpit 

A pilot was needed for the aircraft for a number of reasons. Some 
of these reasons were; performance of landing maneuvers, and 
handling any problems that might occur during the high speed 
ionization communications blackout. The cockpit was designed to 
accommodate both the pilot and the various computers needed to 
record the experimental flight data. The cockpit was sized to an 
average length, width and height of 11’, 5’, and 5’ respectively. 
Included in this space were provisions for cockpit instrumentation, 
an oxygen system, an air conditioning system, 200 lbs. of on-board 
flight computers, emergency equipment, and an ejection seat. A 
glass window will be included for use at velocities up to Mach 
three. For higher speeds a metal covering plate will protect the 
observation window. The pilot will be wearing a pressure suit, so 
ejection may be accomplished for speeds up to Mach two (at high 
altitude). Weights of the various cockpit furnishings are shown in 
Appendix H . 

7 . 2 Landing Gear 

Due to an inadequate allotment of time in the previous quarter, one 
of the areas that needed to be addressed this quarter was landing 
gear. Reference 10 was used for this analysis. This analysis 
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encompassed many aircraft variables such as landing speed, center 
of gravity location, and length of aircraft. Possible locations 
for the landing gear were taken from the drawing of the aircraft 
which indicated where there was adequate room to store the gear. 
The location for the nose gear was set at 17 feet from the nose of 
the aircraft. The location of the main gear (skids) was set at 65 
feet from the nose. A decision to use two tires for the nose gear 
was made to aid in the controllability of the aircraft upon 
landing. Skids were used for the rear because there was not 
adequate space in the aircraft to house conventional wheeled gear. 
Drawings of the nose gear and the main gear during operation are 
shown in Figures 7-1 & 7-2 respectively. Reference 10 used an 
energy method to determine the sizing of the shock absorber. The 
total amount of energy that needed to be absorbed at touch-down was 
determined. The energy absorbed by each set of gear was set to the 
total value that needed to be absorbed. The reason for this was 
because one set will touch down before the other one does. 

There were many variables in determining the size of our shock 
absorbers. Table 7—1 shows some of the values that were used in 
the analysis. 
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Table 7-1: 


Landing Gear Design Factors 


Vertical landing speed = 10 ft/sec 

Landing load factor = 8 

Absorbing medium efficiency = 0.85 (liquid springs) 


After these values in Table 7-1 were decided upon, the analysis 
could be completed. Values of the sizes of the shock absorbers are 
shown in Table 7-2. The length of the struts shown in the drawings 
is arbitrary as long as the shock absorber stroke is met, and there 
is no danger of the aircraft contacting the landing surface upon 
full compression of the shock absorber. 

Table 7-2: Landing Gear Parameters 


Nose gear (2 tires) 

Required shock absorber length (stroke) 
Required shock absorber diameter 
Energy absorbed 


.992 feet 
.311 feet 
66770 BTU 


Main gear (skids) 

Required shock absorber length (stroke) = .737 feet 
Required shock absorber diameter = .354 feet 
Energy absorbed = 66770 BTU 
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8 . 0 Environmental Aspects 


The major reason that supersonic or hypersonic aircraft cannot fly 
over land is because of a sonic boom which creates overpressures. 
This sonic boom is due to the sudden rise and fall of pressure 
waves. The maximum increase in pressure is called the 
overpressure, and is measured in terms of psf . An overpressure of 
less then 1.0 psf is considered relatively safe. 


There are many ways that the overpressure of an aircraft can be 
calculated. A method by Morris was chosen. The rise in pressure 
is due to either volume or lifting effects, whichever is greater. 
It was determined that the rise in pressure due to lift would be 
greater, and the following equation was used: 


dP L - 0 .363 K s 


N p 0 


K, 


(M 2 - 1) 

M 


0.375 


• JR lw 


-0.25 rj-0.5 u-0.75 


Pg 


A plot of the overpressure related to the time in flight can be 
seen in Figure 8-1. 

Because a sonic boom can break glass, frighten animals, and harm 
the human eardrum, it was decided that the flight path of STB-White 
must take place over water. Even though the overpressures created 
are almost always less then 1.0 psf, it was decided to keep the 
flight far from land for safety purposes. 
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9 . 0 Economic Analysis 


An economic analysis was performed using the methods outlined in 
Reference 9. The calculations were done using 1974 dollars, and 
then converted using a consumer price index value determined from 
Fundamentals of Macroeconomics by L. Fleisher. Three basic 
parameters were used for estimating the aircraft cost: the AMPR 
weight, the maximum speed in knots and the quantity of aircraft 
produced. The aircraft had an AMPR weight of 20100 lbs, and a 
maximum speed of 5935 knots. Calculations were performed to 
determine the cost of producing either one or two aircraft, and 
these are compared in Table 9-1. The cost of producing an 
additional proto-type is about 1.15 billion dollars in addition to 
the base cost of 5.45 billion dollars. From this information, it 
would appear that a second aircraft would indeed be produced to aid 
in the SCRAMjet testing. 

A pie chart of the aircraft cost breakdown is shown in Figure 9-1. 
From this figure it can be seen that the majority of the cost for 
the aircraft is in tooling costs. A pie chart of the entire 
development cost breakdown is shown in Figure 9-2. In this figure 
it is apparent that the main expenditure is for the engineering 
required to develop the aircraft. Once the engineering and tooling 
has been completed for one aircraft, it is relatively easy to 
produce other aircraft at a much reduced cost. This is the reason 
for only an approximate increase of 20 % for an additional aircraft. 
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This cost estimation may be somewhat optimistic considering that 
many new technologies will need to. be developed to complete the 
aircraft. A listing of economic analysis is shown in Appendix J. 


Table 9-1: Aircraft Development Costs 


TOTAL DEVELOPMENT COST TO PRODUCE ONE AIRCRAFT 


AIRFRAME ENGINEERING 


$2,631,030,667 

DEVELOPMENT SUPPORT 


$1,958,839,157 

AIRCRAFT 


$805,571 ,066 

ENGINES & AVIONICS 

$17,202,500 


MANUFACTURING LABOR 

$126,512,434 


MATERIAL & EQUIPMENT 

$16 , 298,046 


TOOLING 

$629,111,470 


QUALITY CONTROL 

$16,446,616 


FLIGHT TEST OPERATIONS 


$55,340,978 

TEST FACILITIES 


$0 


TOTAL COST = 

$5,450,781,868 


TOTAL DEVELOPMENT COST TO PRODUCE TWO AIRCRAFT 


AIRFRAME ENGINEERING 
DEVELOPMENT SUPPORT 
AIRCRAFT 

ENGINES & AVIONICS 
MANUFACTURING LABOR 
MATERIAL $ EQUIPMENT 
TOOLING 

QUALITY CONTROL 
FLIGHT TEST OPERATIONS 
TEST FACILITIES 


$2,986,856,785 

$2,489,746,643 

$996,027,024 

$17 , 202 , 500 
$181,916, 851 
$28,219,628 
$745 , 038 , 854 
$23 , 649 , 191 

$134,482 , 648 
$0 


TOTAL COST 


$6 , 607 ,113,101 


Aircraft Cost 


Aircraft C c- s t 
Figure 9-1 
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Total Development Cost 


Total Development Cost 

Figure 9-2 
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Development Support 
36 % 


CONCI.USIONS 


With design work complete, the configuration of the STB-White has 
been finalized. Implementing propulsion, waverider studies with 
aerodynamics, cockpit, and landing systems design, the aircraft 
configuration is complete. The aircraft utilizes two General 
Electric turbof anramjet engines and five full-size SCRAMjet 
combustion cells, also from General Electric, occupying a Mach 10.0 
modified waverider configuration. The turbof anram jets use 6,460 
pounds of liquid hydrogen to carry the aircraft from Mach 0.8 at an 
altitude of 40,000 feet to Mach 6.0 at 80,000 feet. The five 
SCRAMjets are then activated accelerating the aircraft from Mach 
6.0 to Mach 10.0 at 100,000 feet. The SCRAMjets hold the aircraft 
at this altitude and Mach number for two minutes during which 
important data is collected and after which all engines are shut 
down for a glide descent and landing. The total fuel expenditure 
for the SCRAMjet phase is 7428 pounds of liquid hydrogen. Since 
the modified Mach 10.0 waverider has roughly 6,500 cubic feet 
internal volume, the 13,828 total pounds of fuel at 3,121 cubic 
feet of fuel volume can be contained in an airf rarae/fuel tank 
integrated structure. Problems involving the pressurization of the 
fuel and building the airframe from materials resistant to liquid 
hydrogen corrosion are solved. The aircraft is configured with a 
diamond-shape planform with a length of 85 feet, maximum width of 
48 feet, and height (including vertical tails) of 14 feet. A 
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planform area of 2,100 square feet gives the aircraft an aspect 
ratio of roughly 1.1. The cockpit will include avionics, pilot 
ejection systems, a pilot pressure suit, and a cockpit cooling 
system. The landing system will include a front nose wheel for 
directional control upon landing and rear sleds. Furthermore, the 
landing gear configurations will be as simple and compact as 
possible to conserve space and reduce weight. With these design 
aspects complete, the STB-White design is complete. 
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SCRAMjet Information 
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Appendix A: 
SCRAMjet Information 


The following material does not represent the present conditions of 
the aircraft. It has been included to show how the process of 
selection was reached and should not be taken as an accurate 
depiction of fuel weights, accelerations, drags or any thing else. 
It was not revised because the employment of five engines had been 
finalized and only the corresponding data was revised throughout 
the quarter. A graph of total fuel versus number of SCRAMjets can 
be seen in Figure A-l . This figure should put the importance of 
the iteration procedure in perspective. The following page 
contains a section of the spreadsheet which is discussed in the 
preceding sections. This page contains all of the information for 
the cases of 3 to 6 engines. It is broken up in to the four 
divisions previously mentioned. Under the column titled MACH # , 

the first line in each division represents the range over which the 
data is used. Halfway down the column in each division is the Mach 
number that was used to calculate data in that division. The 
accuracy of this spread sheet could be greatly increased if data 
was available allowing 40 or 50 divisions as opposed to the 4 
divisions the author has developed. The acceleration is dependent 
upon the mass of the airplane. The amount of fuel used in the 
previous section by five engines is subtracted from the total mass 
of the airplane. This allows for a constantly changing mass. The 
mass of the fuel from five engines was subtracted because that 
appeared to be the number of engines that would be used. 
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3/6/90 

PJG 

SCRAMJET 

INFORMATION 



# ENG 

NET F 

ISP 

TIME 

WEIGHTf 

ALT (ft) 

MACH # ACC 

5 

8500 

3000 

0.023904 

1219.112 

80000 

6-6.75 8.436466 

6 

8500 

3000 

0.014498 

887.2916 

80000 

13.90974 

3 

8500 

3000 

-0.08034 

-2458.47 

80000 

-2.51008 

4 

8500 

3000 

0.068057 

2776.737 

80000 

2 . 963189 

5 

11500 

1800 

0.011144 

1281.648 

80000 

18.09519 

6 

11500 

1800 

0.007908 

1091 . 363 

80000 

25.50021 

3 

11500 

1800 

0.061387 

4235.731 

80000 

6 3.285146 

4 

11500 

1800 

0.018864 

1735.550 

80000 

10.69016 

5 

12500 

1250 

0.009461 

1703.044 

80000 

21 .31476 

6 

12500 

1250 

0.006867 

1483.464 

80000 

29.36370 

3 

12500 

1250 

0.038656 

4174.899 

80000 

5.216891 

4 

12500 

1250 

0.015201 

2189.082 

80000 

13.26582 

5 

6600 

2497 

0.065029 

3093.899 

88875 

6.75-8.256.300590 

6 

6600 

2497 

0.038641 

2206.160 

88875 

10.60306 

3 

6600 

2497 

-0.17780 

-5075.60 

88875 

-2.30436 

4 

6600 

2497 

0.205054 

7804.718 

88875 

1 .998113 

5 

9337 

1570 

0.026916 

2881 .417 

88875 

15.22171 

6 

9337 

1570 

0.019228 

2470.016 

88875 

7.5 21.30841 

3 

9337 

1570 

0.134409 

8633.005 

88875 

3.048308 

4 

9337 

1570 

0.044851 

3841.056 

88875 

9.135010 

5 

10236 

1106 

0.022571 

3760.223 

88875 

18.15195 

6 

10236 

1106 

0.016504 

3299.394 

88875 

24.82471 

3 

10236 

1106 

0.085244 

8520.465 

88875 

4.806457 

' 4 

10236 

1106 

0.035692 

4756.803 

88875 

11.47920 

5 

5300 

2150 

0.055838 

2477.693 

97000 

8.25-9.254.974620 

6 

5300 

2150 

0.032531 

1732.166 

97000 

8.538843 

3 

5300 

2150 

-0.12896 

-3433.58 

97000 

-2.15382 

4 

5300 

2150 

0.196950 

6991.270 

97000 

1.410396 

5 

8000 

1350 

0.019766 

2108.374 

97000 

14.05330 

6 

8000 

1350 

0.014293 

1829.623 

97000 

19.43326 

3 

8000 

1350 

0.084344 

5398.029 

97000 

9 3.293382 

4 

8000 

1350 

0.032026 

2732.937 

97000 

8.673342 

5 

8500 

980 

0.017654 

2756.187 

97000 

15.73453 

6 

8500 

980 

0 .012949 

2426 . 060 

97000 

21 .45074 

3 

8500 

980 

0.064567 

6048.268 

97000 

4.302125 

4 

8500 

980 

0.027726 

3463.039 

97000 

10.01833 

5 

4673 

1850 

0.050613 

2301.251 

101340 

9.25-10 4.116151 

6 

4673 

1850 

0.028438 

1551.642 

101340 

7.325629 

3 

4673 

1850 

-0.09046 

-2468.02 

101340 

-2.30280 

4 

4673 

1850 

0.229777 

8357.850 

101340 

0.906673 

5 

7086 

1247 

0.016797 

1718.125 

101340 

12.40255 

6 

7086 

1247 

0.012063 

1480.717 

101340 

10 17.26931 

3 

7086 

1247 

0.078055 

4790.301 

101340 

2.669035 

4 

7086 

1247 

0.027645 

2262.178 

101340 

7.535793 

5 

8025 

908 

0.013331 

2120.854 

101340 

15.62714 

6 

8025 

908 

0.009855 

1881 .442 

101340 

21 . 13881 

3 

8025 

908 

0.045252 

4319.429 

101340 

4.603788 

4 

8025 

908 

0.020595 

2621.166 

101340 

10.11546 


WEIGHT FOR 2 MINUTE CRUISE 
595.6027 
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NUMBER OF ENGINES VS. FUEL WEIGHT 
FOR SCRAMJET OPERATION 
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NUMBER OF ENGINES 


Appendix B: 

Turbof anramjet Inlet Code 
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c ********************************************************************* 
program TURBOFAN 

********************************************************************* 
AERONAUTICAL & ASTRONAUTICAL ENGINEERING * 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

C 

C 

C 

c 

c 

c 

c 

C 

C 

c 

C 

c 


416 HYPERSONIC DESIGN 
SPRING QUARTER 1990 

AUTHORS: GOUHIN, P. & WHITE, H. 


PURPOSE: 


to AID IN THE DESIGN OF AN INLET FOR A TURBOFANRAM JET 
ENGINE ON A HYPERSONIC TEST BED VEHICLE. THIS 
CODE ALLOWS FOR THE ADDITION OF SEVERAL RAMPS, ON 
THE UNDERSIDE OF THE VEHICLE, to INDUCE SHOCKS 
to INCREASE PRESSURE RATIOS AND REDUCE FLOW 
VELOCITY. 

PRESENTLY: THE PROGRAM IS DESIGNED TO CONSIDER THREE RAMPS ON 
THE UNDER SIDE OF THE VEHICLE. FROM THIS POINT THE 
FLOW CURVES INTO THE COMBUSTION CHAMBER AT SUBSONIC 
SPEEDS. 


INPUT: 


AS DEFINED IN THE programfset up for five 
INTERACTIVE VARIABLES FOR FREESTREAM MACH 
AND RAMP ANGLES ARE INPUT FROM KEYBOARD. 


ramps ) 
NUMBER 


OUTPUT: DATA IS WRITTEN to TF.DAT 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


c ********** ***************************************************** ****** 
implicit real *8 (a-h,l-z) 


integer i , j , k 

open( 1 , f ile= ’TF. DAT’ , STATUS= ’ NEW ’ ) 

30 format ( 5X, ’VALUES FROM PROGRAM TURBOFAN FOR’) 

40 format (3X, ’FREESTREAM MACH # = ’ , F6 . 3 ) 

50 format (3x, ’betal=’ , f 6 . 3 , 3x , ’ M2 = ’ , f 5 . 3 , 3x , ’Pt2/Ptl=’ ,f6.4 

F , 3X , ’P2/Pl=’ , F6 . 4 , 3X , ’ R02/R01 = ’ , F7 . 4 ) 

60 f ormat ( 3x , ’ be t a 2 =’ , f 6 . 3 , 3x , ’ M3 = ’ ,f5.3,3x, ’Pt3/Pt2= ,f6.4 

F, 3X , ’P3/P2=’ , F6 . 4 , 3X , ’ R03/R02= ’ , F7 . 4 ) 

7 0 format (3x, ’beta3=’ , f 6 . 3 , 3x , ’ M4 = ’ ,f5.3,3x, ’Pt4/Pt3= ,f6.4 

F,3X, ’ P4 /P3 = ’ , F6 . 4 , 3X , ’R04/R03=’ ,F7.4) 

80 format (3x, ’beta4=’ , f 6 . 3 , 3x , ’ M5 = ’ , f 5 . 3 , 3x , ’Pt5/Pt4=’ ,f6.4 

F,3X, ’ P5 /P4 = ’ ,F6.4,3X, ’R05/R04=’ ,F7.4) 

90 format(3x, ’ M6 = ’ ,f5.3,3x, ’Pt6/Pt5=’ ,f6.4 

F,3X, ’ P6/P5= ’ , F8.4 ,3X, ’R06/R05=’ , F7.4 ) 

100 format (3x, ’beta6=’ , f 6 . 3 , 3x , ’ M7= ’ , f 5 . 3 , 3x , ’Pt7/Pt6=’ ,f6.4 
F , 3 X , ’ P7 /P6= ’ , F6.4 ,3X, ’R07/R06=’ ,F7.4 ) 

Cl 05 formaU 1 7X , ’ M8= ’ , F5 . 3 , 1 7X , ’ P8/P7= ’ F6 . 4 , ’R08/R07= ’ , F7 .4 ) 

110 format ( 5X , ’ THE TOTAL PRESSURE RECOVERY IS ’ , FG . 4 ) 

120 format! 5X , ’THE STATIC PRESSURE RECOVERY IS ' , F8 . 4 ) 

130 formaU 5X, ’THE DENSITY RATIO IS ’,F8.4) 

140 format ( 5X , ’THE ITERATED BETA 2 IS ’ , F6 . 3 ) 

150 format! 5x, ’THE MASS FLOW RATE IS ’.F10.4) 

160 format! 5x , ’THE PRESSURE IN THE COMBUSTION CHAMBER IS ’ , FI 0 . 4 ) 
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C OTHERWISE KNOWN AS PI . 

pi = 3.14159265358979323846264338327950288 

C THE FREESTREAM FLOW IS DEFINED 

WRITE (*,*)" ENTER FREESTREAM MACH NUMBER!!" 

READ( 5 , * )M1 

WRITE (*,*)" ENTER AMBIENT TEMPERATURE!!" 

READ ( 5 , * )T1 

WRITE (*,*)" ENTER FREESTREAM DENSITY!!!!” 

READ ( 5 , * ) ROl 
R01=R01*32 . 2 

WRITE (*,*)" ENTER AMBIENT PRESSURE !!!!!!!!!" 

READ( 5 , * ) PI 
AREA=3059 . /144 . 0 
write (1,30) 
write ( 1 , 4 0 ) M 1 

THETA1 IS THE HALF-WEDGE NOSE ANGLE OF THE VEHICLE 
WRITE (*,*)" ENTER WEDGE ANGLE OF VEHICLE!!!!" 

READ(5,*)THETA1 
thetal = THETAl*pi/180 . 

BETA IS NEVER EXPECTED to EXCEED 70 DEG. AND THUS IT IS 
ITERATED DOWN FROM THIS INITIAL guess 
betal = 70.*pi/180. 

LEFT AND RIGHT REFER to THE L.H.S & R.H.S. OF THE THETA- 
BETA-MACH NUMBER EQUATION 
leftl = dtan(thetal) 
rt 1 = 1.0 

5 betal=betal*dabs(leftl/rtl)**0.1 

rta=( (ml**2*dsin( betal ) **2 ) -1 ) / ( ml**2* ( 1 . 4+dcos ( 2*betal ) ) +2 ) 
rtb=2/dtan( betal ) 
rtl=rta*rtb 

if THE TWO SIDES ARE NOT EQUAL WITHIN THIS RANGE then THE 
ITERATION PROCESS MUST continue AND A NEW BETA IS CHOSEN 
if ( dabs ( le f t 1 -r t 1 ) . g t . 1 d- 6 ) goto 5 
C THE NORMAL FREESTREAM MACH NUMBER IS CALCULATED 
mnl =ml *ds in ( be tal ) 

C THE NORMAL MACH NUMBER AFTER THE SHOCK IS CALCULATED 
mn2=dsqrt ( ( 1 + . 2*mnl ** 2 . )/(l . 4*mnl **2 . - . 2 ) ) 

C THE PRESSURE RATIO ACCROSS THE SHOCK IS CALCULATED 

pt2ptl=(mnl/mn2)*( ( 2+ ( 0 . 4*mn2**2 ) )/(2+( .4*mnl**2) ) )**3.00 
p2pl=l . 0+2 .8/2 . 4* (mnl** 2 . 0-1 .0 ) 
c FOR A CALORICALLY PERFECT GAS. 

ro2rol =2 . 4*mnl**2 . /( 0 . 4 *mn 1 * * 2 . + 2 . ) 


T2T1=P2P1/R02R01 

TTA=( (1 + 1 . 4 *mn 1 * * 2 . )/( 1 + 1 . 4*mn2**2. ) ) **2* ( mn2/mn 1 )**2 
TTB= ( 1 . + . 2*mn2**2 ) / ( 1 .+.2*mnl**2) 


TT2TT1 =TTA*TTB 

THE MACH NUMBER BEHIND THE SHOCK IS CALCULATED 


m2 = mn2 /d s i n ( be t a 1 - 1 he t a 1 ) 
betal = betal * 1 80 . /pi 

write) 1 , 50 ) betal ,m2,pt2ptl,P2Pl, R02 RO 1 * * # ^ 

C VALUES ARE DEFINED AND THE FIRST RAMP IS ANALYZED 


******** 
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WRITE (*,*)" ENTER FIRST RAMP ANGLE!!!!" 

READ ( 5 , * ) THETA2 
theta2 = THETA2*pi/180 . 
beta2 = 70*pi/180. 

Ieft2 = dtan(theta2) 
rt2 = 1.0 

6 beta2=beta2*dabs (Ieft2/rt2)**0. 1 

rta=( (m2**2*dsin(beta2)**2)-l)/(m2**2*(1.4+dcos(2*beta2) ) + 2) 

rtb=2/dtan( beta2 ) 

rt2=rta*rtb 

if ( dabs ( lef t2-rt2 ) . gt . Id-5 ) goto 6 
mn2=m2*dsin( beta2 ) 

mn3=dsqrt ( ( 1+ . 2*mn2**2 . ) / ( 1 . 4*mn2**2 . - . 2 ) ) 

pt3pt2 = ( mn2/mn3 ) * ( ( 2+ ( 0 . 4*mn3 **2 ) ) / ( 2+ ( . 4*mn2**2 ) ) ) ** 3 . 00 

p3p2=l .0+2.8/2.4*(mn2**2.-l .0) 

ro3ro2=2 . 4*mn2**2 . / ( 0 . 4*mn2**2 . +2 . ) 

T3T2=P3P2/R03R02 

TTA= ( (1+1 . 4*mn2**2 . )/( 1+1 . 4*mn3**2 . ) ) **2* ( mn3/mn2 ) **2 
TTB= ( 1 . +. 2*mn3**2 )/( 1 . + . 2*mn2**2 ) 

TT3TT2 =TTA*TTB 

m3 = mn3/dsin( beta2-theta2 ) 

beta2 = beta2*180 . /pi 

write(l,60) beta2 ,m3 ,pt3pt2 ,P3P2 , R03R02 
c ********************************************************************* 
C VALUES ARE DEFINED AND THE SECOND RAMP IS ANALYZED 
WRITE (*,*)" ENTER SECOND RAMP ANGLE!!!!" 

READ( 5 , * )THETA3 
theta3 = THETA3*pi/180 . 
beta3 = 70.*pi/180. 

Ieft3 = dtan(theta3) 
rt3 =1.0 

7 beta3=beta3*dabs(left3/rt3)**0 . 1 

rta=( (m3**2*dsin(beta3 ) ** 2 )-l )/ ( m3**2* ( 1 . 4+dcos ( 2 *beta3 ) ) + 2 ) 

rtb=2/dtan( beta3 ) 

rt3=rta*rtb 

if ( dabs ( lef t3-rt3 ). gt . Id-4 ) goto 7 
mn3=m3*dsin( beta3 ) 

mn4 =dsqrt ( ( 1+ . 2*mn3**2 . )/( 1 . 4*mn3**2 . - . 2 ) ) 

Pt4pt3=(mn3/mn4 )*( ( 2+ ( 0 . 4*mn4**2 ) )/(2+( .4*mn3**2) ) )**3.00 
p4p3=l .0+2. 8/2. 4*(mn 3**2. 0-1.0) 
ro4ro3=2.4*mn3**2./(0.4*mn3**2.+2.) 

T4T3=P4P3/R04R03 

TTA= ( (1 + 1 . 4*MN3**2 . )/( 1 + 1 .4*MN4**2. ) ) **2* ( MN4 /MN3 )**2 
TTB= ( 1 . + . 2*MN4**2 )/( 1 . +.2*MN3**2) 

TT4TT3=TTA*TTB 

m4 = mn4 /ds i n ( be ta3 - the ta3 ) 

beta3 = beta3* 1 80 . /pi 

write! 1 ,70) be ta3 ,m4 ,pt4pt3,P4P3, R04 R03 

C********************************************************************* 

C VALUES ARE DEFINED AND THE THIRD RAMP IS ANALYZED 
WRITE (*,*)" ENTER THIRD RAMP ANGLE!!!!" 

READ( 5 , * )THETA4 
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theta4 = THETA4*pi/180 . 

C theta4 = 3.0*pi/180. 

beta4 = 70.*pi/180. 

Ieft4 = dtan(theta4) 
rt4 =1.0 

8 beta4=beta4*dabs ( lef t4/rt4 ) **0 . 1 .. . 

rta=( (m4**2*dsin(beta4 )**2)-l )/(m4**2*( 1 . 4+dcos ( 2*beta4 ) )+2 ) 

rtb=2/dtan( beta4 ) 
rt4 =rta*rtb 

WRITE (*,*) LEFT4 , RT4 
if (dabs( Ieft4-rt4) .gt.ld-4) goto 8 
mn4=m4*dsin(beta4 ) 

mn5=dsqrt ( ( 1 + . 2*mn4**2 . ) / ( 1 . 4*mn4**2 . - . 2 ) ) 

pt5pt4= ( mn4/mn5 ) * ( ( 2+ ( 0 . 4*mn5**2 ) ) / ( 2+ ( . 4*mn4**2 ) ) ) **3 . 00 
p5p4=1.0+2.8/2.4*( mn4 ♦*2.0— 1.0) 
ro5ro4 = 2.4 *inn4 * * 2 . / ( 0 . 4 *inn4 * * 2 . + 2 . ) 

T5T4=P5P4/R05R04 , 

i PTA= ((1 + 1.4 *MN4**2 . )/( 1 + 1 .4*MN5**2. ) ) **2* ( MN5/MN4 ) ** 2 

TTB= (l.+.2*MN5**2)/(l.+. 2*MN4**2 ) 

TT5TT4 =TTA*TTB 

m5 = mn5/dsin( beta4-theta4 ) 
beta4 = beta4*180 . /pi 

write(l,80) beta4 , ra5 , pt5pt4 , P5P4 , R05R04 0 . 0 .******** 

C****************** *************************************************** 

C VALUES ARE DEFINED AND THE FOURTH RAMP IS ANALYZED 
C theta5 = 3.0*pi/180. 

C beta5 = 70.*pi/180. 

C left5 = dtan(theta5) 

C rt5 = 1.0 

C 9 beta5=beta5*dabs(left5/rt5)**0.2 

C rta=( (m5**2*dsin(beta5)**2)-l )/( m5**2* ( 1 . 4+dcos ( 2*beta5 ) )+2) 

G rtb=2/dtan ( beta5 ) 

C rt5=rta*rtb 

C if ( dabs ( lef t5-rt5 ). gt . Id-4 ) goto 9 

C mn5=m5*dsin( beta5 ) 

C ran6=dsqrt( ( 1+ . 2*mn5**2 . )/( 1 . 4*mn5**2 .-. 2 ) ) 

C P t6pt5=(mn5/mn6)*((2+(0.4*mn6**2))/(2+(.4*mn5**2)))**3.00 

C p6p5=l. 0+2. 8/2. 4*(mn 5**2. 0-1 .0) 

C r o6ro5 = 2.4*mn5**2./(0.4*mn5**2.+2. ) 

C T6T5=P6P5/R06R05 

C TTA= ( ( 1 + 1 . 4*MN5**2 . )/( 1 + 1 . 4*MN6**2 . 

C TTB= (l.+.2*MN6**2)/(l.+.2*MN5**2) 

C TT6TT5 = TTA*TTB 

C m6 = mn6/dsin( beta5-theta5 ) 

C beta5 = beta5* 1 80 . /p i 

C write! 1 ,90) beta5,n6,pL6pt5,P6P5,R06RO5 

c * ***************************** *************************************** 


) **2* ( MN6/MN5 )**2 


C 

c 

c 

c 

c 


VALUES ARE DEFINED AND THE FIFTH RAMP IS ANALYZED 
theta6 = 3.0*pi/180. 
beta6 = 70.*pi/180. 

Ieft6 = dtan(theta6) 
rt6 = 1.0 
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C 10 beta6=beta6*dabs ( le f t6/rt6 ) **0 . 2 

C rta= ( (m6**2*dsin(beta6)**2)-l ) / ( m6**2* ( 1 . 4+dcos ( 2*beta6 ) )+2 ) 

C rtb=2/dtan( beta6 ) 

C rt6=rta*rtb 

C if ( dabs ( lef t6-rt6 ) . gt . Id-4 ) goto 10 

C mn6=m6*dsin( beta6 ) 

C mn7=dsqrt ( ( 1+ . 2*ran6**2 . ) / ( 1 . 4*mn6**2 . - . 2 ) ) 

C Pt7pt6=(mn6/mn7)*( ( 2+ ( 0 . 4*mn7**2 ) ) / ( 2+ ( .4*mn6**2) ) )**3.00 

C p7p6=1.0+2.8/2.4*( mn6**2 .0-1 .0) 

C ro7ro.6 = 2 . 4*mn6**2 . / ( 0 . 4*mn6**2 . + 2 . ) 

C T7T6=P7P6/R07R06 

C TTA= ( ( 1+1 .4*MN6**2 . ) / ( 1 + 1 . 4*MN7**2 . ) ) **2* ( MN7/MN6 ) **2 

C TTB= ( 1 . + . 2*MN7**2 )/( 1 .+ . 2*MN6**2 ) 

C TT7TT6=TTA*TTB 

C m7 = mn7/dsin( beta6-theta6 ) 

C beta6 = beta6* 1 80 . /pi 

C write(l,100) beta6 , m7 , pt 7pt6 , P7P6 , R07R06 

c ******************************************************************* 

C NORMAL SHOCK OCCURS AND VALUES ARE CALCULATED FOR COMBUSTION 

m6=dsqrt( ( 1+ . 2*m5**2 . ) / ( 1 . 4*m5**2 . - . 2 ) ) 
pt6pt 5 = ( m5 / m6 ) * ( ( 2+ ( 0 . 4 *m6**2 ) ) / ( 2 + ( . 4*m5**2 )))**3.00 
p6p5=l .0+2. 8/2. 4*(m5**2. 0-1.0) 
ro6ro5=2.4*m5**2./(0.4*ra5**2.+2. ) 

T6T5=P6P5/R06R05 

TTA= ( ( 1 + 1 . 4*M5**2 . )/( 1 + 1 .4*M6**2 . ) ) **2* ( M6/M5 ) **2 
TTB= ( 1 . + . 2*M6**2 )/( 1 . + . 2*M5**2 ) 

TT6TT5=TTA*TTB 

write (1,90) m6 ,pt6pt5,P6P5, R06R05 
C THE FOLLOWING PART OF THE PROGRAM WAS WRITTEN FOR THE 

C EXPANSION WAVE THAT WOULD BE GENERATED AS THE FLOW IS 

C TURNED BACK INTO THE SCRAMJET COMBUSTION CHAMBER. 

C HOWEVER, A COWL WAS ADDED AT THE BOTTOM OF THE SCRAMJET 

C MODULE TO GET AWAY FROM THIS EXPANSION. THEREFORE, AT THIS 

C POINT IN TIME THE EXPANSION WAVE IS NOT VALID. 

c ******************************************************************* 
C THE CONDITIONS FOR THE EXPANSION WAVE ARE CALCULATED 

C USING THE PRANDTL- MEYER function. 

C PART = DSQRT (0.4/2.4*(M7**2.-l .0) ) 

C PARTI =DATAN( PART) 

C NU1=SQRT( 6.0) 

C NU2 = DATAN ( DSQRT ( M 7 ** 2 . - 1 . ) ) 

C NU7 = PART 1 *NU1 -NU2 

C THETA7=21 . *PI/180 . 

C NU8=THETA7+NU7 

C NU8=76 . 92 

C GNU8=NU8-5 . 0 

C M8 = 3 . 0 

C 1 2 M8=M8* ( NU8/GNU8 ) ** . 3 

C NU3 =SQRT ( 6. 0 ) *DATAN ( DSQRT ( 0.4/2. 4* (M8** 2.-1 .0) ) ) 

C NU4 = DATAN ( DSQRT ( M8**2 . - 1 . ) ) 

C GNU 8= ( NU3-NU4 ) * 1 80 . /PI 

C write(*,* )GNU8,NU8,M8 
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non 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


if ( DABS ( GNU8-NU8 ) . gt . ID-5 ) goto 12 

NOTE THAT THE TOTAL PRESSURES AND TEMPERATURES do NOT 
CHANGE THROUGH A PRANDTL -MEYER EXPANSION WAVE. 

T8T7=l/( ( 1+ . 4/2 . 0*M8**2 . ) / ( 1+ . 4/2 . *M7**2 . )) 

P8P7= 1 / ( ( ( 1+ .4/2 . *M8**2. ) / ( 1 + . 4/2 . *M7**2 . ) )**3 .5 ) 

R08R07=P8P7/T8T7 

write (1,105) M8 , P8P7 , R08R07 

THE FOLLOWING VALUES ARE SET TO ONE SO THOSE RAMPS COULD 
BE "C"ED OUT AND THE TOTAL RATIOS ARE STILL CALCULATED 
AT THE END OF THE PROGRAM. FOR EACH ADDITIONAL RAMP 
THAT IS ANALYZED, THE CORRESPONDING NUMBERS MUST BE DELETED 
SO THAT THE VALUES ARE NOT RESET TO ONE. 

P8P7 = 1 . 0 
R08R07 = 1 . 0 
T8T7=1 . 0 
PT7PT6 =1.0 
P7P6 = 1 . 0 
TT7TT6=1 . 0 
R07R06=1 . 0 
T7T6 =1.0 
PT6PT5=1 .0 
P6P5 = 1 . 0 
TT6TT5 = 1 . 0 
R06R05 = 1 . 0 
T6T5=1 .0 
PT5PT4=1 . 0 
P5P4=1 . 0 
TT5TT4=1 . 0 
R05R04 = 1 . 0 
T5T4=1 . 0 

PTRATIO=PT2PTl *PT3PT2 *PT4PT3 *PT5PT4 *PT6PT5*PT7PT6 
PRATIO= 1 . /(P2P1*P3P2*P4P3*P5P4*P6P5*P7P6*P8P7) 

TTRATIO=TT2TTl *TT3TT2*TT4TT3*TT5TT4*TT6TT5*TT7TT6 

R0RATI0=R02R01 *R03R02*R04R03*R05R04*R06R05*R07R06*R08R07 

TRATIO=T2Tl*T3T2*T4T3*T5T4*T6T5*T7T6*T8T7 

R05=R01 *RORATIO 

T5=Tl*TRATIO 

THE SPEED OF SOUND IS CALCULATED 
A=DSQRT( 1. 4*1716. *T5) 

THE VELOCITY RIGHT BEFORE THE COMBUSTION CHAMBER IS CALCULATED 


V=M5*A 

C MASS FLOW RATE ENTERING THE COMBUSTION CHAMBER. 

MDOT=R05*V*AREA 
PRESSURE=P1 /PRAT IO 
write! 1,110) PT RATIO 
write! 1 , 120 )PRATIO 
wr i te ( 1 , 1 30 ) RORATI O 
wri te( 1 , 1 50 )MDOT 
WRITE! 1 , 160 )PRESSURE 

Jr****************************************************************** 

POSITIONING THE SHOCK WAVE SO THAT IT IMPINGES ON THE LIP OF 
THE SCRAMJET. THIS IS AN ITERATIVE PROCESS. 
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o o n 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 


11 


THIS IS THE CALCULATION OF THE DOWNWARD ANGLE OF THE 
HORIZONTAL to EACH SHOCK WAVE. 

TH1 = THETA 1 
TH2=TH1 +THETA2 
TH3=TH2+THETA3 
TH4=TH3+THETA4 
TH5 =TH4 + THETA5 
TH6=TH5+THETA6 
BT1 =BETA1 
BT2=TH1+BETA2 
BT2=BT2*PI/180 . 

BT3=TH2+BETA3 

BT4=TH3+BETA4 

BT5=TH4+BETA5 

BT6=TH5+BETA6 

USING SIMULAR TRIANGLES, THE RAMPS CAN BE PLACE SO THAT 
THE OBLIQUE SHOCKS HIT AT THE LIP OF THE ENGINE. 

THESE ARE THE LENGTHS OF THE ORIGINAL TRIANGLE THAT 
DEFINES THE AREA BEFORE THE FIRST RAMP BEGINS. 

LI = 4 7 . 0 
HI =8 . 287 

DEFINITION OF THE DIMENSIONS OF THE FIRST RAMP. 

L2=50 . 2-L1 
H2=10 . 6-H1 
H2L2 = 1 . 0 
AH2L2 = 1 . 0 

BT2=BT2*DABS(H2L2/AH2L2 )**0 . 3 
H2L2 = H2/L2 
AH2L2=DTAN ( BT2 ) 

if ( DABS (H2L2-AH2L2) .gt. ID- 6) goto 11 
EBETA2=BT2-TH1 
EBETA2=EBETA2*1 80 . /PI 
WRITE ( 1 , 140 )EBETA2 
stop 
end 


109 


Appendix C: 


SCRAMjet Inlet Program 
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Q********************************************************************* 


program INLET 

c ***************************** ************************************* 
c AERONAUTICAL & ASTRONAUTICAL ENGINEERING 

c 416 HYPERSONIC DESIGN 

c SPRING QUARTER 1990 


c 

c 

c 

c 

C 

c 

C 

c 

c 


AUTHORS: GOUHIN, P. & WHITE, H. 

PURPOSE: to AID IN THE DESIGN OF AN INLET FOR A SCRAMJET 
ENGINE ON A HYPERSONIC TEST BED VEHICLE. THIS 
CODE ALLOWS FOR THE ADDITION OF SEVERAL RAMPS, ON 
THE UNDERSIDE OF THE VEHICLE, to INDUCE SHOCKS 
to INCREASE PRESSURE RATIOS AND REDUCE FLOW 
VELOCITY. 


*** 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

C 

c 

C 


PRESENTLY: THE PROGRAM IS DESIGNED TO CONSIDER ONE RAMP ON THE * 
UNDER SIDE OF THE VEHICLE. THE SECOND RAMP TAKES IN- * 
TO ACCOUNT THE COWL ANGLE. THE THIRD RAMP TAKES INTO * 
ACCOUNT THE SECOND ANGLE ON THE COWL. FROM THIS POINT* 
THE FLOW GOES DIRECTLY INTO THE COMBUSTION CHAMBER * 
FOR A CLEARER EXPLANATION WITH FIGURES SEE THE FINAL * 

REPORT FOR THE WHITE GROUP ( AAE 416 SPRING 1990) * 

* 

INPUT: AS DEFINED IN THE program(set up for five ramps) * 

INTERACTIVE VARIABLES FOR FREESTREAM MACH NUMBER * 

AND RAMP ANGLES ARE INPUT FROM KEYBOARD. * 

* 

OUTPUT: DATA IS WRITTEN to INLET . DAT * 


c 


* 


c ********************************************************************* 


implicit real *8 (a-h,l-z) 
integer i , j , k 

open( 1 , file= * INLET . DAT ’ , STATUS= ’ NEW ’ ) 

30 ’ format) 5X, ’VALUES FROM PROGRAM INLET FOR’) 

40 format (3X, ’FREESTREAM MACH # = ’ , F6 . 3 ) 

50 f or mat ( 3x, ’ be t al = ’ , f 6 . 3 , 3x , ’M2=’ ,f5.3,3x, ’Pt2/Ptl = ’ ,f6.4 

F,3X, ’P2/P1=’ , F6 . 4 , 3X , ’ R02 /RO 1 = ’ , F7 . 4 ) 

60 format (3x, ’beta2=’ , f 6 . 3 , 3x , ’ M3= ’ , f 5 . 3 , 3x , ’Pt3/Pt2=’ ,f6.4 

F,3X, ’ P3/P2 = ’ , F6 . 4 , 3X , ’R03/R02=’ ,F7.4) 

70 f ormat ( 3x , ’ be ta3 = ’ , f 6 . 3 , 3x , ’ M4 = ’ , f 5 . 3 , 3x , ’Pt4/Pt3=* ,f6.4 

F,3X, ’ P4 /P3 = ’ , F6 . 4 , 3X , ’ R04 /R03 = ’ , F7 . 4 ) 

80 format(3x, ’beta4=’ , f 6 . 3 , 3x , ’M5=’ , f 5 . 3 , 3x , ’Pt5/Pt4=’ ,f6.4 

F,3X, ’ P5 /P4 = ’ ,F6.4,3X, ’ R05 /R04 = ’ , F7 . 4 ) 

90 format(3x, ’beta5=' , f 6 . 3 , 3x , ’M6=’ , f 5 . 3 , 3x , ’Pt6/Pt5=’ ,f6.4 

F,3X, ' P6/P5 = ’ , F6 . 4 , 3X , ’R06/R05=’ , F7.4 ) 

100 format ( 3x, ’ b«*ta6= ’ , f 6 . 3 , 3\ , ’ M7= ’ , f 5 . 3 , 3x , ’Pt7/Pt6= ’ , f6.4 
F,3X, ’ P7/P6 = ' ,F6.4,3X, ’R07/R06= ’ ,F7.4 ) 

Cl 05 format ( 1 7X , ’M8= ’ , F5. 3 , 1 7X , ’P8/P7= ’ F6 . 4 , ’R08/R07 = ’ , F7 . 4 ) 
110 format ( 5X , 'THE TOTAL PRESSURE RECOVERY IS ’,F6.4) 

120 formatf 5X , ’THE STATIC PRESSURE RECOVERY IS ’ , F8 . 4 ) 

130 format) 5X, ’THE DENSITY RATIO IS ’ , F8 . 4 ) 
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140 format ( 5X, ’THE ITERATED BETA2 IS ’ , F6 . 3 ) 

150 format ( 5x, ’THE MASS FLOW RATE IS ’ , F8 . 4 ) 

160 format (5x, ’THE PRESSURE IN THE COMBUSTION CHAMBER IS ’ ,F10.4) 
C OTHERWISE KNOWN AS PI. 

pi = 3.14159265358979323846264338327950288 
C THE FREESTREAM FLOW IS DEFINED 

WRITE (*,*)" ENTER FREESTREAM MACH NUMBER!!" 

READ( 5 , * ) Ml 

WRITE (*,*)" ENTER AMBIENT TEMPERATURE!!" 

READ( 5 , * )T1 

WRITE (*,*)" ENTER FREESTREAM DENSITY!!!!" 

READ ( 5 , * ) ROl 
ROl =R01 * 3 2 . 2 

WRITE (*,*)" ENTER AMBIENT PRESSURE!!!!!!!!!" 

READ ( 5 , * )P1 

AREA = 8 . 0*30 . 0/144 . 0 

write (1,30) 

write ( 1 , 40 ) M 1 

C THETA 1 IS THE HALF-WEDGE NOSE ANGLE OF THE VEHICLE 
WRITE (*,*)" ENTER WEDGE ANGLE OF VEHICLE!!!!" 

READ( 5 , * )THETA1 
thetal = THETA 1 *pi / 180. 

BETA IS NEVER EXPECTED to EXCEED 70 DEG. AND THUS IT IS 
ITERATED DOWN FROM THIS INITIAL guess 
betal = 70.*pi/180. 

LEFT AND RIGHT REFER to THE L.H.S & R.H.S. OF THE THETA- 
BETA-MACH NUMBER EQUATION 
leftl = dtan(thetal) 
rt 1 = 1.0 

5 betal = betal *dabs (leftl/rtl)**0.3 

rta= ( (ml**2*dsin( betal )**2)-l)/(ml**2*(l . 4+dcos ( 2* betal ) ) +2 ) 
rtb=2/dtan ( betal ) 
rtl =rta*rtb 

if THE TWO SIDES ARE NOT EQUAL WITHIN THIS RANGE then THE 
ITERATION PROCESS MUST continue AND A NEW BETA IS CHOSEN 
if ( dabs ( le f t 1 -rt 1 ) . g t . 1 d-6 ) goto 5 
C THE NORMAL FREESTREAM MACH NUMBER IS CALCULATED 
mnl =ml *ds in ( betal ) 

C THE NORMAL MACH NUMBER AFTER THE SHOCK IS CALCULATED 
mn2=dsqrt ((l+.2*mnl**2.)/(1.4*mnl**2.-.2)) 

C THE PRESSURE RATIO ACCROSS THE SHOCK IS CALCULATED 

pt2ptl=(mnl/mn2 )*( ( 2+ ( 0 . 4*mn2**2 ) )/(2+( .4*mnl**2) ) )**3.00 
p2pl = 1 .0+2. 8/2. 4*(mn 1**2. 0-1.0) 
c FOR A CALORICALLY PERFECT GAS. 

ro2 ro 1 = 2 . 4 *mn 1 * * 2 . / ( 0 . 4 *mn 1 * * 2 . + 2 . ) 

T2T1=P2P1 /RO2R01 

TTA= ( (1 + 1. 4*mn 1**2. ) / ( 1 ♦ 1 . 4 *mn2 * * 2 . ) ) * *2 * ( mr>2 /mn 1 )**2 
TTB= ( 1 .+ .2*mn2**2 )/( 1 ,+.2*mnl**2 ) 

TT2TT1 =TTA*TTB 

C THE MACH NUMBER BEHIND THE SHOCK IS CALCULATED 
m2 = mn2/dsin{ betal -thetal ) 
betal = betal * 1 80 . /pi 
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write (1,50) betal i m2 , pt2pt 1 , P2P1 , R02R01 

C VALUES ARE DEFINED AND THE FIRST RAMP IS ANALYZED 
WRITE (*,*)" ENTER FIRST RAMP ANGLE!!!!" 

READ ( 5 , * )THETA2 
theta2 = THETA2*pi/180 . 
beta2 = 70*pi/180. 

Ieft2 = dtan(theta2) 
rt2 = 1.0 

6 beta2=beta2*dabs(left2/rt2)**0.3 

rta= ( (m2**2*dsin(beta2 ) **2 ) -1 ) / ( m2**2* ( 1 . 4 + dcos ( 2*beta2 ) )+2 ) 

rtb=2/dtan(beta2) 

rt2=rta*rtb 

if ( dabs ( lef t2-rt2 ) . gt . Id-5 ) goto 6 
mn2=m2*dsin( beta2 ) 

mn3=dsqrt ( ( 1 + . 2 *mn2** 2 . ) / ( 1 . 4 *mn2 **2 . - . 2 ) ) 

pt3pt2= ( mn2/mn3 )*( (2+(0.4*mn3**2) ) / ( 2+ ( . 4*mn2**2 ) ) )**3.00 

p3p2=l .0+2. 8/2. 4*(mn 2**2. -1 .0) 

ro3ro2=2 . 4*mn2**2 . /( 0 . 4*mn2**2 . +2 . ) 

T3T2=P3P2/R03R02 

TTA= ((1+I.4*mn2**2. ) / ( 1+1 . 4*mn3**2 . ) ) **2* ( mn3/mn2 ) **2 
TTB= ( 1 . + . 2*mn3**2 )/(!.+. 2*mn2**2 ) 

TT3TT2 =TTA*TTB 

m3 = mn3/dsin( beta2-theta2 ) 

beta2 = beta2*180 . /pi 

write( 1 , 60 ) beta2 , m3 , pt3pt2 , P3P2 , R03R02 

c ************************************************************* ******** 
C VALUES ARE DEFINED AND THE SECOND RAMP IS ANALYZED 
WRITE (*,*)" ENTER COWL ANGLE!!!!" 

READ( 5 , * )THETA3 
theta3 = THETA3*pi/180 . 
beta3 = 70.*pi/180. 

Ieft3 = dtan(theta3) 
rt3 = 1.0 

7 beta3=beta3*dabs (Ieft3/rt3)**0 . 3 

rta= ( ( m3**2*dsin(beta3 )**2 )-l ) / ( m3** 2* ( 1 . 4+dcos ( 2*beta3 ) ) +2 ) 

rtb=2/dtan(beta3) 

rt3=rta*rtb 

if ( dabs ( lef t3-rt3 ). gt . Id-4 ) goto 7 
mn3=m3*dsin(beta3 ) 

mn4 =dsqrt ( ( l+.2*mn3**2. )/( 1 . 4 *mn3**2 . - . 2 ) ) 

P t4pt3=(mn3/mn4 )*( ( 2+ ( 0 . 4*mn4**2 ) )/(2+( .4*mn3**2) ) )**3.00 
p4 p3= 1 .0+2. 8/2. 4*(mn3**2. 0-1.0) 
ro4 ro3 = 2 . 4*mn3**2 . / ( 0 .4*mn3**2.+2. ) 

T4T3=P4P3/R04R03 

TTA= ( (1 + 1 . 4*MN3**2 . ) / ( 1 + 1 . 4 *MN4 * * 2 . > ) **2* ( MN4/MN3 )**2 
TTB= ( 1 . + . 2*MN4**2 ) /( 1 . + . 2*MN3**2 ) 

TT4TT3 =TTA*TTB 

m4 = mn4 /ds i n ( be La3- t he La3 ) 

beta3 = be ta3 * 1 80 . /p i 

write! 1 ,70) beta3,m4,pt4pt3,P4P3, R04 R03 

Qt ************************************************************* ******* 
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C VALUES ARE DEFINED AND THE THIRD RAMP IS ANALYZED 
WRITE (*,*)" ENTER SECOND COWL ANGLE!!!!" 

READ ( 5 , * ) THETA4 
theta4 = THETA4*pi/180 . 

C ■ theta4 = 3.0*pi/180. 

beta4 = 70.*pi/180. 

Ieft4 = dtan(theta4) 
rt4 = 1.0 

8 beta4=beta4*dabs ( Ieft4/rt4)**0.1 

rta= ( (m4**2*dsin(beta4 )**2)-l )/(m4**2*( 1 .4+dcos( 2*beta4 ) )+2 ) 

rtb=2/dtan ( beta4 ) 

rt4=rta*rtb 

WRITE ( * , * ) LEFT4 , RT4 
if (dabs(left4-rt4) .gt.ld-4) goto 8 
mn4=m4*dsin( beta4 ) 

mn5=dsqrt ( ( 1+ . 2*mn4**2 . ) / ( 1 . 4*ron4**2 . - . 2 ) ) 

Pt5pt4=(mn4/mn5 )*( ( 2+ ( 0 . 4*mn5**2 ) )/(2+( .4*mn4**2 ) ) )**3 .00 
p5p4=l .0+2.8/2.4*( mn4**2 .0-1 .0) 
ro5ro4=2 . 4*mn4**2 . /(0 . 4*mn4**2 . +2 . ) 

T5T4=P5P4/R05R04 

TTA= ( ( 1 + 1 . 4*MN4**2 . ) / ( 1 + 1 . 4 *MN5**2 . ) ) **2* ( MN5/MN4 ) ** 2 
TTB= ( 1 . + . 2*MN5**2 )/(!.+. 2*MN4**2 ) 


TT5TT4 =TTA*TTB 

m5 = mn5/dsin( beta4-theta4 ) 

beta4 = beta4* 1 80 . /pi 

write(l,80) beta4 , m5 , pt5pt4 , P5P4 , R05R04 

********************************************* *********** 

C VALUES ARE DEFINED AND THE FOURTH RAMP IS ANALYZED 
C theta5 = 3 . 0*pi/180 . 

C beta5 = 70.*pi/180. 

C left5 = dtan(theta5) 

C rt 5 = 1 . 0 

C 9 beta5=beta5*dabs(left5/rt5)**0.2 

C rta= ( (m5**2*dsin(beta5)**2)-l)/(m5**2*(1.4+dcos(2*beta5) )+2) 

C rtb=2/dtan( beta5 ) 

C rt5=rta*rtb 

C if ( dabs ( lef t5-rt5 ) . gt . Id-4 ) goto 9 

C mn5=m5*dsin( beta5 ) 

C mn6 = dsqrt ( (l+.2*mn5**2. ) / ( 1 . 4 *mn5 **2 . - . 2 ) ) 

C P t6pt5=(mn5/mn6 )*( ( 2+( 0 . 4*mn6**2 ) )/( 2+( .4*mn5**2 ) ) )**3 . 00 

C P 6p5=l .0+2. 8/2. 4*(mn 5**2. 0-1 .0) 

C ro6ro5=2.4*mn5**2./(0.4*mn5**2.+2. ) 

C T6T5=P6P5/R06R05 

C TTA=( (1 + 1 . 4 *MN5**2 . )/( 1 + 1 . 4*MN6**2. ) )**2*(MN6/MN5 )**2 

C TTB= ( 1 . + . 2*MN6**2 )/( 1 . + . 2*MN5**2 ) 

C TT6TT5=TTA*TTB 

C m6 = mn6 /ds i n ( be La5 - t he ta5 ) 

C be ta5 = be ta5* 1 80 . /pi 

C write(l,90) beta5 , m6 , pt6pt5 , P6P5 , R06R05 ************ 

C VALUES ARE DEFINED AND THE FIFTH RAMP IS ANALYZED 
C theta6 = 3.0*pi/180. 


114 


C beta6 = 70.*pi/180. 

C left6 = dtan(theta6) 

C rt6 = 1.0 

C 10 beta6=beta6*dabs ( le f t6/rt6 ) **0 . 2 

C rta= ( (m6**2*dsin(beta6)**2)-l)/(m6**2*(l. 4+dcos ( 2*beta6 ) )+2 ) 

C rtb=2/dtan( beta6 ) 

C rt6=rta*rtb 

C if ( dabs ( lef t6-rt6 ) . gt . Id-4 ) goto 10 

C mn6=m6*ds in ( beta6 ) 

C mn7 =dsqrt((l+. 2*mn6** 2 . ) / ( 1 . 4 *mn6**2 . - . 2 ) ) 

C pt7pt6=(mn6/mn7)*( ( 2+ ( 0 . 4*mn7**2 ) )/(2+( . 4*mn6**2 ) ) )**3.00 

C p7p6= 1 .0+2.8/2.4*( mn6**2 .0-1.0) 

C ro7ro6 = 2 . 4*mn6**2 . /( 0 . 4*mn6**2 . +2 . ) 

C T7T6=P7P6/R07R06 

C TTA= ( (1+1 .4*MN6**2. ) / ( 1 + 1 . 4 *MN7**2 . ) ) **2* ( MN7/MN6 ) **2 

C TTB= ( 1 . + . 2*MN7**2 )/(l.+.2*MN6**2) 

C TT7TT6=TTA*TTB 

C ra7 = ran7/dsin( beta6-theta6 ) 

C beta6 = beta6*180 . /pi 

C write(l,100) beta6 , m7 , pt7pt6 , P7P6 , R07R06 

c ******************************************************************* 
C THE FOLLOWING PART OF THE PROGRAM WAS WRITTEN FOR THE 

C EXPANSION WAVE THAT WOULD BE GENERATED AS THE FLOW IS 

C TURNED BACK INTO THE SCRAM JET COMBUSTION CHAMBER. 

C HOWEVER, A COWL WAS ADDED AT THE BOTTOM OF THE SCRAM JET 

C MODULE TO GET AWAY FROM THIS EXPANSION. THEREFORE, AT THIS 

C POINT IN TIME THE EXPANSION WAVE IS NOT VALID. 

************* ************************************************** 
C THE CONDITIONS FOR THE EXPANSION WAVE ARE CALCULATED 

C USING THE PRANDTL- MEYER function. 

C PART=DSQRT (0.4/2.4*(M7**2.-l .0) ) 

C PART 1=DATAN( PART) 

C NU1 =SQRT (6.0) 

C NU2=DATAN ( DSQRT( M7**2 . - 1 . ) ) 

C NU7=PART1 *NU1 -NU2 

C THETA7=21 . *PI/180 . 

C NU8 = THETA7 + NU7 

C NU8=76 . 92 

C GNU8=NU8-5 . 0 

C M8=3 . 0 

C 1 2 M8=M8* ( NU8/GNU8 )** . 3 

C NU3=SQRT( 6 . 0 ) * DAT AN ( DSQRT ( 0.4/2.4*(M8**2.-1.0))) 

C NU4=DATAN( DSQRT( M8**2 . -1 . ) ) 

C GNU8 = ( NU3 -NU4 ) * 1 80 . /PI 

C write!*, * )GNU8 , NU8 , M8 

C if ( DABS(GNU8-NU8 ) .gt . lD-f. ) goto 12 

C NOTE THAT THE TOTAL PRESSURES AND TEMPERATURES do NOT 

C CHANGE THROUGH A PRANDTL- MEY ER EXPANSION WAVE. 

C T8T7 = 1 / ( (l + .4/2.0*M8**2. )/(!♦. 4/2. *M7»* 2. ) ) 

C P8P7 = 1 / ( ( ( l + .4/2.*M8**2. )/( l + .4/2.*M7**2. ) )**3.5 ) 

C R08R07 = P8P7/T8T7 

C write(l,105) M8 , P8P7 , R08R07 
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THE FOLLOWING VALUES ARE SET TO ONE SO THOSE RAMPS COULD 
BE "C"ED OUT AND THE TOTAL RATIOS ARE STILL CALCULATED 
AT THE END OF THE PROGRAM. FOR EACH ADDITIONAL RAMP 
THAT IS ANALYZED, THE CORRESPONDING NUMBERS MUST BE DELETED 
SO THAT THE VALUES ARE NOT RESET TO ONE. 

P8P7= 1 . 0 
R08R07= 1 . 0 
T8T7= 1 .0 
PT7PT6 =1.0 
P7P6= 1 . 0 
TT7TT6=1 . 0 
R07R06= 1 . 0 
T7T6 =1.0 
PT6PT5=1 .0 
P6P5=1 . 0 
TT6TT5 = 1 . 0 
R06R05 = 1 . 0 
T6T5= 1 .0 
C PT5PT4 = 1.0 

C P5P4=1 . 0 

C TT5TT4=1 . 0 

C R05R04 = 1.0 

C T5T4 = 1 . 0 

PTRATIO=PT2PTl*PT3PT2*PT4PT3*PT5PT4*PT6PT5*PT7PT6 

PRAT I 0=1 . / ( P2P1*P3P2*P4P3*P5P4*P6P5*P7P6*P8P7 ) 

TTRATIO=TT2TTl*TT3TT2*TT4TT3*TT5TT4*TT6TT5*TT7TT6 

R0RATI0=R02R01 *R03R02*R04R03*R05R04*R06R05*R07R06*R08R07 

TRATIO=T2Tl*T3T2*T4T3*T5T4*T6T5*T7T6*T8T7 

R05=R01 *RORATIO 

T5=T1 *TRATIO 

C THE SPEED OF SOUND IS CALCULATED 

A=DSQRT ( 1. 4*1716. *T5) 

C THE VELOCITY RIGHT BEFORE THE COMBUSTION CHAMBER IS CALCULATED 

V=M5*A 

C MASS FLOW RATE ENTERING THE COMBUSTION CHAMBER. 

MDOT=R05*V*AREA 
PRESS URE=P1 /PRAT IO 
write ( 1 , 110 )PTRATIO 
write ( 1 , 120 )PRATIO 
wri te ( 1 , 130) RORATIO 
wr i te ( 1 , 150 )MDOT 
WRITE ( 1 , 160 )PRESSURE 

C********** ************** ******************************************* 

C POSITIONING THE SHOCK WAVE SO THAT IT IMPINGES ON THE LIP OF 

C THE SCRAMJET. THIS IS AN ITERATIVE PROCESS. 

C 

C THIS IS THE CALCULATION OF THE DOWNWARD ANGLE OF THE 

C HORIZONTAL to EACH SHOCK WAVE. 

TH 1 =THETA1 

C TH2=TH 1 +THETA2 

C TH3=TH2+THETA3 

C TH4=TH3+THETA4 
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TH5=TH4+THETA5 
TH6=TH5+THETA6 
BT1 =BETA1 
BT2=TH1 +BETA2 
BT2=BT2*PI/180 . 

C BT3=TH2+BETA3 

C BT4=TH3+BETA4 

C BT5=TH4+BETA5 

C BT6=TH5+BETA6 

C USING SIMULAR TRIANGLES, THE RAMPS CAN BE PLACE SO THAT 

C THE OBLIQUE SHOCKS HIT AT THE LIP OF THE ENGINE. 

C THESE ARE THE LENGTHS OF THE ORIGINAL TRIANGLE THAT 

C DEFINES THE AREA BEFORE THE FIRST RAMP BEGINS. 

Ll=47 . 0 
Hl=8 . 287 

C DEFINITION OF THE DIMENSIONS OF THE FIRST RAMP. 

L2=50 . 2-L1 
H2= 1 0 . 6-H1 
H2L2=1 . 0 
AH2L2=1 . 0 

11 BT2=BT2*DABS(H2L2/AH2L2 )**0 . 3 
H2L2 = H2/L2 
AH2L2=DTAN( BT2 ) 

if ( DABS ( H2L2-AH2L2 ) . gt . ID-6 ) goto 11 
EBETA2=BT2-TH1 
EBETA2=EBETA2* 1 80 . /PI 
WRITE ( 1 , 140 )EBETA2 
stop 
end 
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Appendix D: 
Inlet Code Output 



™™™SchT^,’ 1' ,!so ™ for 

betal= 3 ? 2.000 

beta2=36.*686 M 3 -i* 7 nf ^/Ptl-O. 9993 
beta3=43 .197 SJ-J ’ S&E* 2 " 0 * 99 ®* 

beta4=52 . 804 KS -1 206 Pt4/Pt3«0.9956 

^-6 = 0.8 38 pfr Pt5/Pt4=0. 9957 

THE TOTAL P pRESSnpr'f| 22 P6 /P5= 1.5301 

THE STATIC p£ff^ r ^£ 0VERY TS 0.98il 1 

s ss^^nir “ ° ■ 2065 

THE treses „ the ioMBuIJioS’Jffiuan, Is 


P2/P1-1.2069 
P3/P2=l . 3022 
P4/P3=l .4205 
P5/P4=l .4174 


' — • A / *4 Kfj 

R06/R05= 1.3520 


4913 . 0079 


R02/R01= 
R03/R02= 
P04/'R03 = 
P05/R04 = 


F0E 

S: BE bbs*® 

beta4=4 8 : ill • 3 02 ? PtVPt3=0. 8972 

M6=0 .766 Pt6/Pt5=n 3 Q7T, Pt5 / p b4=0 . 9197 

THE TOTAL PRESSURE ’rECOVFRY^tc 5 ^ 1 * 9303 
THE STATIC PRrco^- „£ 0VERY Is 0.6844 
THE DeSI Tv^ffc^i^^RV I£ 0 .01« 

the mass eeow rate IS 


P2/P1=1 .6199 
P3/P2=3 .1307 
P4/P3=2 . 8003 
P5/P4=2 . 5678 


' ~ w w / o KU 

R06/R05= 1.5865 


R02/P01= 
R03/R02= 
R04/R03= 
R05/R04 = 


1.143 
1.2 06 
1.2:3 
1.22: 


.406 
. 166 
.022 
.914 


«-U <.u 0» Oi 


VALUES FROM PROGRAM INLET FOR 
FREESTREAM MACH # « 6.000 


betal 

beta2 

beta3 

beta4 

THE 

THE 

THE 

THE 

THE 


12.372 

15.349 

17.452 

16.676 


M2=5. 449 
M3=4 . 696 
M4=4 . 058 
M5=3 .795 


Pt2/Pt 1=0 .9818 
Pt3/Pt2=0 .9470 
Pt4/Pt3=0. 9561 
Pt5/Pt4=0 .9957 


TOTAL PRESSURE RECOVERY IS 0.8851 
STATIC PRESSURE RECOVERY IS 0.0826 
DENSITY RATIO IS 5.7335 
MASS FLOW RATE IS 142.6178 
PRESSURE IN THE COMBUSTION CHAMBER IS 


P2/P1=1 .7614 
P3/P2=2 .2608 
P4/P3=2 . 1478 
P5/P4=l .4155 


703 . 6571 


R02/R01= 
R03/R02= 
R04/R03- 
R05/R04 = 


VALUES FROM PROGRAM INLET FOR 
FREESTREAM MACH # = 7.000 

betal=ll . 018 M2=6.281 Pt2/Ptl=0 . 9723 

beta2=12 .739 M3=5.541 Pt3/Pt2=0. 9619 

beta3=13 . 947 M4=4.942 Pt4/Pt 3=0 . 9727 

beta4=14 .443 M5=4.538 Pt5/Pt4=0 . 9893 

THE TOTAL PRESSURE RECOVERY IS 0.8999 
THE STATIC PRESSURE RECOVERY IS 0.0818 
THE DENSITY RATIO IS 5.8035 
THE MASS FLOW RATE IS 135.2251 
_THE PRESSURE IN THE COMBUSTION CHAMBER IS 


P2/P1=1 . 9215 
P3/P2=2 . 0712 
P4/P3=l . 9141 
P5/P4=l . 6057 


560.5341 


R02/R01= 

R03/R02= 

R04/R03= 

R05/R04= 


VALUES FROM PROGRAM INLET FOR 
FREESTREAM MACH - = 8.000 

betal=10 . 019 M2=7.087 Pt2/Ptl=0 . 9602 

beta2=ll .728 M3=6.173 Pt3/Pt2=0 . 9475 

K e * a x = ^’ 723 ^4=5 .313 Pt4/Pt3=0.9406 

beta4=12 . 875 M5=4.970 Pt5/Pt4=0 . 9942 

THE TOTAL PRESSURE RECOVERY IS 0.8509 
THE STATIC PRESSURE RECOVERY IS 0.0618 
THE DENSITY RATIO IS 6.9769 
THE MASS FLOW RATE IS 127.6593 
THE PRESSURE IN THE COMBUSTION CHAMBER IS 


P2/P1=2 . 0935 
P3/P2=2 .2543 
P4/P3=2 . 3355 
P5/P4=l . 4687 


511.9135 


R02/R01= 

R03/R02= 

R04/R03= 

R05/R04= 


VALUES FROM PROGRAM INLET FOR 


FREESTREAM MACH # = 
betal= 9.255 M2=7. 

beta2=10 . 136 M3=6.981 

beta3=ll . 036 M4=6.271 


9.000 

868 Pt2/Ptl=0 .9456 
Pt3/Pt2=0.9619 
Pt4/Pt3=0 .9726 
beta4=ll .963 M5=5.674 Pt5/Pt4=0 . 9794 

THE TOTAL PRESSURE RECOVERY IS 0.8664 
THE STATIC PRESSURE RECOVERY IS 0.0613 
THE DENSITY RATIO IS 7.0507 
THE MASS FLOW RATE IS 120.9392 
THE PRESSURE IN THE COMBUSTION CHAMBER IS 


P2/P1=2 .2776 
P3/P2=2 . 0701 
P4/P3=l . 9168 
P5/P4=l . 8047 


430.4089 


R02/R01= 

R03/R02= 

R04/R03= 

R05/R04= 


VALUES FROM PROGRAM INLET FOR 
FREESTREAM MACH # = 10.000 
betal= 8.653 M2=8.623 Pt2/Ptl=0 . 9283 

beta2= 9.119 M3-7.704 Pt3/Pt2=0.9661 

beta3«ll . 536 M4-6.535 Pt4/Pt3-0. 9164 

beta4=10. 118 M5-6.206 Pt5/Pt4-0. 9968 

THE TOTAL PRESSURE RECOVERY IS 0.8193 
THE STATIC PRESSURE RECOVERY IS 0.0563 
THE DENSITY RATIO IS 7.3747 
THE MASS FLOW RATE IS 114.5754 
THE PRESSURE IN THE COMBUSTION CHAMBER IS 


P2/P1=2 .4742 
P3/P2=2 . 0119 
P4/P3=2 . 6023 
P5/P4-1.3711 


383.7050 


R02/R01= 

R03/R02= 

R04/R03= 

R05/R04= 



1.5816 
1 . 6636 
1.5775 
1.3982 


1 . 6755 
1.7598 
1.8011 
1.3137 


.7717 

.6630 

1.5790 

1.5155 


1 .8698 
1.6315 
1.9313 
1.2517 
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Appendix E: 


SCRAMjet Inlet Representative Configurations 





SCRRMjet Inlets 

(not draLim to scale) 






Appendix F: 
Aerodynamic Data 
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AR = 1 . 1 


alpha =4 
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AR= 1 . 1 


alpha=4 
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AR = 1 . 1 


alpha=2 
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AR= 1 . 1 


alpha=2 



40K 

50K 

60K 

MACH 

L/D 

L/D 

L/D 

0 . 2 

7.111 

6 . 944 

6.768 

0.4 

7.347 

7.190 

7.026 

0.6 

7 .494 

7.343 

7.184 

0 . 8 

4.432 

4.381 

4.326 

1 .0 

1 . 208 

1 .201 

1 .194 

1.2 

2.623 

2.619 

2.615 

1 . 4 

2.910 

2.906 

2.901 

1 . 6 

3.739 

3.729 

3.719 

1.8 

4 .169 

4 .156 

4.141 

2.0 

4 . 648 

4.629 

4.608 

2 . 2 

4 . 987 

4 . 964 

4 . 938 

2.4 

5.299 

5.269 

5.236 

2 . 6 

5.547 

5.514 

5.479 

2 . 8 

5 . 608 

5.567 

5.523 

3 . 0 

5.763 

5.719 

5 . 671 

3.2 

6.027 

5 . 982 

5 . 933 

3.4 

6 . 137 

6.089 

6.037 

3.6 

6. 197 

6.146 

6.091 

3.8 

6 .400 

6.348 

6.291 

4.0 

6.520 

6.466 

6.407 

4 . 2 

6.618 

6.562 

6.502 

4 . 4 

6.725 

6 . 667 

6.606 

4 . 6 

6.885 

6.826 

6 .763 

4 . 8 

6.999 

6.940 

6.875 

5.0 

7.100 

7.039 

6.974 

5 . 2 

7.141 

7.079 

7.013 

5 . 4 

7.257 

7.195 

7 .127 

5.6 

7.311 

7.248 

7.180 

5 . 8 

7.405 

7.341 

7.272 

6.0 

7.427 

7 . 362 

7.293 

6.4 

7.552 

7.486 

7.416 

7.0 

7.893 

7 . 826 

7.754 

7.4 

8.007 

7 . 939 

7 . 866 

8.0 

8.054 

7.985 

7.912 

8 . 4 

8.284 

8.215 

8.141 

9.0 

8.303 

8.234 

8.160 

9 . 4 

8.488 

8.419 

8 . 344 

10.0 

8 . 523 

8.454 

8.379 


70K 

80K 

90K 

100K 

L/D 

L/D 

L/D 

L/D 

6.585 

6 . 395 

6.194 

5 . 993 

6.853 

6.673 

6.482 

6 . 289 

7.016 

6.841 

6.655 

6.466 

4.267 

4.204 

4.135 

4.065 

1 . 187 

1 .178 

1 .169 

1 . 160 

2.611 

2.606 

2.601 

2.595 

2.896 

2.890 

2.884 

2.877 

3.707 

3.694 

3.680 

3.664 

4.125 

4.108 

4.089 

4.069 

4 . 585 

4.561 

4.533 

4 . 505 

4.911 

4.881 

4.848 

4 . 813 

5.201 

5.162 

5.120 

5.075 

5.440 

5.399 

5 . 353 

5 . 304 

5.475 

5.423 

5 . 366 

5 . 307 

5.620 

5.565 

5 . 504 

5.440 

5.880 

5.823 

5 . 760 

5.695 

5.981 

5.920 

5.854 

5.784 

6.032 

5.968 

5.898 

5.824 

6.230 

6.165 

6.092 

6.017 

6.344 

6.276 

6.202 

6.124 

6.437 

6.367 

6.290 

6.210 

6.539 

6.467 

6.389 

6 . 307 

6.696 

6.622 

6.542 

6.459 

6.806 

6.732 

6.650 

6.565 

6.904 

6.828 

6.745 

6.658 

6.941 

6.864 

6.779 

6 . 691 

7.055 

6.977 

6.891 

6.801 

7.106 

7.027 

6.940 

6.850 

7.198 

7.118 

7.030 

6 . 938 

7.217 

7.137 

7.048 

6.956 

7.339 

7.257 

7.167 

7.073 

7 .675 

7 . 591 

7.499 

7 .403 

7 . 787 

7 . 702 

7.609 

7.512 

7.832 

7 . 747 

7.653 

7.555 

8.060 

7.974 

7.879 

7 . 780 

8.079 

7.993 

7.898 

7 .799 

8.263 

8.175 

8.080 

7 . 980 

8.298 

8.211 

8.116 

8.016 


129 



AR = 1 . 1 

t / c = l /12 

alpha =0 


40 K 

50 K 

60 K 

MACH 

CD 

CD 

CD 

0.2 

0.016 

0.017 

0.018 

0.4 

0.015 

0.016 

0.017 

0.6 

0.015 

0.015 

0.016 

0.8 

0.023 

0.023 

0.024 

1 . 0 

0.044 

0.045 

0.045 

1.2 

0.162 

0.163 

0. 163 

1 . 4 

0.135 

0.136 

0.136 

1 . 6 

0.077 

0.078 

0.078 

1.8 

0.061 

0.061 

0.062 

2.0 

0.046 

0.047 

0.047 

2.2 

0.039 

0.039 

0.039 

2.4 

0.031 

0.032 

0.032 

2 . 6 

0.028 

0.029 

0.029 

2 . 8 

0.023 

0.023 

0.024 

3 . 0 

0.021 

0.021 

0.021 

3 . 2 

0.020 

0.020 

0.020 

3 . 4 

0.018 

0.018 

0.019 

3 . 6 

0.017 

0.017 

0.017 

3.8 

0.016 

0.016 

0.016 

4 . 0 

0.015 

0.015 

0.015 

4 . 2 

0.014 

0.014 

0.014 

4 . 4 

0.013 

0.013 

0.013 

4 . 6 

0.012 

0.013 

0.013 

4.8 

0.012 

0.012 

0.012 

5.0 

0.011 

0.011 

0.012 

5 . 2 

0.011 

0.011 

0.011 

5 . 4 

0.010 

0.010 

0.010 

5 . 6 

0.010 

0.010 

0.010 

5 . 8 

0.009 

0.010 

0.010 

6 . 0 

0.009 

0.009 

0.009 

6.4 

0.008 

0.008 

0.009 

7 . 0 

0.008 

0.008 

0.008 

7 . 4 

0.007 

0.007 

0.007 

8 . 0 

0.006 

0.007 

0.007 

8 . 4 

0.006 

0.006 

0.006 

9.0 

0.006 

0.006 

0.006 

9.4 

0.005 

0.005 

0.005 

10.0 

0.005 

0.005 

0.005 


70 K 

80 K 

90 K 

100 K 

CD 

CD 

CD 

CD 

0.018 

0.019 

0.020 

0.021 

0.017 

0.018 

0.019 

0.020 

0.017 

0.017 

0.018 

0.019 

0.024 

0.025 

0.026 

0.027 

0.046 

0.046 

0.047 

0.048 

0.164 

0.164 

0.165 

0.165 

0.137 

0.137 

0.138 

0.138 

0.079 

0.079 

0.079 

0.080 

0.062 

0.062 

0.063 

0.063 

0.047 

0.048 

0.048 

0.048 

0.040 

0.040 

0.040 

0.041 

0.032 

0.033 

0.033 

0.033 

0.029 

0.030 

0.030 

0.030 

0.024 

0.024 

0.024 

0.025 

0.022 

0.022 

0.022 

0.022 

0.020 

0.021 

0.021 

0.021 

0.019 

0.019 

0.019 

0.019 

0.017 

0.017 

0.018 

0.018 

0.016 

0.016 

0.017 

0.017 

0.015 

0.015 

0.016 

0.016 

0.014 

0.014 

0.015 

0.015 

0.014 

0.014 

0.014 

0.014 

0.013 

0.013 

0.013 

0.013 

0.012 

0.012 

0.013 

0.013 

0.012 

0.012 

0.012 

0.012 

0.011 

0.011 

0.011 

0.012 

0.011 

0.011 

0.011 

0.011 

0.010 

0.010 

0.010 

0.011 

0.010 

0.010 

0.010 

0.010 

0.009 

0.009 

0.010 

0.010 

0.009 

0.009 

0.009 

0.009 

0.008 

0.008 

0.008 

0.008 

0.007 

0.007 

0.008 

0.008 

0.007 

0.007 

0.007 

0.007 

0.006 

0.006 

0.006 

0.007 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0 . 006 

0.005 

0.005 

0.005 

0.005 
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AR= 1 . 1 


alpha=0 



4 OK 

50K 

60K 

MACH 

L/D 

L/D 

L/D 

0 . 2 

8 . 006 

7.710 

7.410 

0.4 

8.426 

8.139 

7.846 

0.6 

8.675 

8.393 

8.104 

0 . 8 

5 .557 

5.444 

5.326 

1.0 

1.617 

1 . 600 

1.583 

1 . 2 

2.954 

2.947 

2.940 

1 . 4 

3.414 

3.406 

3.397 

1.6 

4.055 

4.038 

4.021 

1.8 

4.429 

4 .408 

4 . 385 

2.0 

4 . 706 

4.679 

4.649 

2 . 2 

4 . 935 

4.903 

4.868 

2.4 

4 . 992 

4.955 

4.915 

2 . 6 

5 . 207 

5.167 

5 . 123 

2.8 

4 .914 

4 .870 

4.823 

3 . 0 

4 . 991 

4 . 945 

4.896 

3.2 

5 .272 

5 .224 

5.171 

3.4 

5 . 289 

5 .239 

5 . 185 

3 . 6 

5.251 

5 . 200 

5 .145 

3 . 8 

5 . 436 

5.383 

5 . 326 

4.0 

5 . 503 

5.449 

5.391 

4 . 2 

5 . 546 

5.491 

5.432 

4 . 4 

5 . 606 

5 .551 

5.491 

4 .6 

5 .739 

5 . 682 

5.621 

4.8 

5.814 

5.757 

5.695 

5.0 

5 .875 

5.817 

5 . 755 

5 . 2 

5 . 865 

5.807 

5.744 

5 . 4 

5 . 950 

5.891 

5 .828 

5 . 6 

5 . 963 

5 . 904 

5.841 

5 . 8 

6.025 

5 . 965 

5 . 902 

6 . 0 

6 . 006 

5.947 

5.883 

6 . 4 

6.068 

6 . 009 

5.945 

7 . 0 

6.335 

6.275 

6.210 

7 . 4 

6 .394 

6 .334 

6 . 269 

8 . 0 

6 . 363 

6 . 304 

6 . 239 

8 . 4 

6.550 

6.490 

6 .425 

9 . 0 

6.506 

6.447 

6 . 383 

9.4 

6 .650 

6.590 

6 . 526 

10.0 

6.633 

6.574 

6.511 


7 OK 

80K 

90K 

100K 

L/D 

L/D 

L/D 

L/D 

7.105 

6.798 

6.483 

6.177 

7.547 

7.244 

6.931 

6.626 

7.808 

7.506 

7.194 

6.889 

5.201 

5.071 

4.932 

4 . 792 

1.564 

1.543 

1 . 521 

1.498 

2.932 

2.924 

2.915 

2 . 904 

3.387 

3.376 

3 . 363 

3 . 351 

4.001 

3.981 

3.957 

3.933 

4.361 

4.334 

4.304 

4.272 

4.617 

4.582 

4 . 544 

4 . 503 

4.831 

4 . 790 

4 . 746 

4.699 

4.871 

4.824 

4 . 772 

4.718 

5 .076 

5.026 

4 . 970 

4.912 

4 . 772 

4.717 

4.657 

4.594 

4.842 

4 . 785 

4 . 722 

4 . 656 

5.115 

5.055 

4.989 

4.920 

5.127 

5.065 

4 . 997 

4 . 926 

5.086 

5.023 

4.953 

4.881 

5.265 

5.200 

5.128 

5.053 

5.329 

5.262 

5 . 189 

5.113 

5.369 

5.301 

5 .226 

5.149 

5 .426 

5 . 358 

5 . 282 

5 . 204 

5.556 

5.486 

5 .409 

5.329 

5.629 

5.558 

5.480 

5 . 399 

5.688 

5.616 

5 . 537 

5.456 

5.677 

5.605 

5.526 

5.444 

5 . 760 

5.687 

5 . 608 

5 . 525 

5.773 

5 . 700 

5.620 

5.537 

5.833 

5 . 760 

5 .679 

5 . 596 

5.815 

5.742 

5 . 662 

5 .578 

5 .877 

5 . 803 

5 . 722 

5 . 639 

6.139 

6.064 

5 . 982 

5 .896 

6.199 

6.123 

6.041 

5 . 955 

6.170 

6.096 

6.015 

5 .930 

6.355 

6.279 

6.197 

6.111 

6 . 314 

6 . 240 

6.159 

6.075 

6.456 

6 . 382 

6 . 300 

6.215 

6.442 

6.369 

6 . 289 

6.205 
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AR = 1 . 1 


alpha = 

-2 


40 K 

50 K 

60 K 

MACH 

CD 

CD 

CD 

0.2 

0.010 

0.011 

0.012 

0.4 

0.009 

0.010 

0.011 

0.6 

0.009 

0.009 

0.010 

0.8 

0.012 

0.012 

0.013 

1.0 

0.020 

0.020 

0.021 

1 . 2 

0.100 

0.100 

0.101 

1.4 

0.077 

0.078 

0.078 

1.6 

0.051 

0.051 

0.051 

1.8 

0.041 

0.042 

0.042 

2.0 

0.033 

0.034 

0.034 

2.2 

0.029 

0.029 

0.029 

2.4 

0.025 

0.025 

0.025 

2 . 6 

0.023 

0.023 

0.023 

2.8 

0.020 

0 . 020 

0 . 020 

3 . 0 

0.018 

0.018 

0.018 

3 . 2 

0.017 

0.017 

0.017 

3 . 4 

0.016 

0.016 

0.016 

3 . 6 

0.015 

0.015 

0.015 

3 . 8 

0.014 

0.014 

0.014 

4.0 

0.013 

0.013 

0.013 

4.2 

0.012 

0.012 

0.013 

4 . 4 

0.012 

0.012 

0.012 

4 . 6 

0.011 

0.011 

0.011 

4.8 

0.011 

0.011 

0.011 

5.0 

0.010 

0.010 

0.010 

5.2 

0.010 

0.010 

0.010 

5.4 

0.009 

0 . 009 

0.009 

5 . 6 

0.009 

0.009 

0.009 

5.8 

0.009 

0 . 009 

0.009 

6.0 

0.008 

0.008 

0.008 

6 . 4 

0.008 

0.008 

0.008 

7 . 0 

0 . 007 

0 . 007 

0.007 

7 . 4 

0.007 

0.007 

0.007 

8.0 

0.006 

0.006 

0.006 

8.4 

0.006 

0.006 

0.006 

9 . 0 

0.005 

0.005 

0.005 

9.4 

0.005 

0.005 

0.005 

10.0 

0 . 005 

0.005 

0 . 005 


70 K 

80 K 

90 K 

100 K 

CD 

CD 

CD 

CD 

0.012 

0.013 

0.014 

0.015 

0.011 

0.012 

0.013 

0.014 

0.011 

0.011 

0.012 

0.013 

0.013 

0.014 

0.015 

0.015 

0.021 

0.022 

0.023 

0.023 

0.101 

0.102 

0.102 

0.103 

0.079 

0.079 

0.080 

0.080 

0.052 

0.052 

0.053 

0.053 

0.042 

0.04 3 

0.043 

0.043 

0.034 

0.035 

0.035 

0 . 036 

0.030 

0.030 

0.030 

0.031 

0.026 

0.026 

0.026 

0.027 

0.023 

0.024 

0.024 

0.024 

0.020 

0.021 

0.021 

0.021 

0.019 

0.019 

0.019 

0 .020 

0.017 

0.018 

0.018 

0.018 

0.016 

0.016 

0.017 

0.017 

0.015 

0.015 

0.016 

0.016 

0.014 

0.014 

0.015 

0.015 

0.013 

0.014 

0.014 

0.014 

0.013 

0.013 

0.013 

0.013 

0.012 

0.012 

0.012 

0.013 

0.011 

0.012 

0.012 

0.012 

0.011 

0.011 

0.011 

0.011 

0.010 

0.011 

0.011 

0.011 

0.010 

0.010 

0.010 

0.010 

0.010 

0.010 

0.010 

0.010 

0.009 

0.009 

0.009 

0.010 

0.009 

0.009 

0.009 

0 . 009 

0.009 

0.009 

0 . 009 

0.009 

0.008 

0.008 

0 . 008 

0 . 008 

0.007 

0.007 

0 . 007 

0.008 

0.007 

0.007 

0.007 

0.007 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0.005 

0.006 

0.006 

0.006 

0.005 

0.005 

0 . 005 

0.005 

0.005 

0.005 

0.005 

0.005 
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7“H 

<r—4 

II 

« 

< 

t/c=l/12 


alpha= 


4 OK 

50K 

60K 

MACH 

L/D 

L/D 

L/D 

0.2 

7 . 793 

7.343 

6.905 

0.4 

8.457 

7.994 

7.540 

0 . 6 

8.846 

8.376 

7.913 

0.8 

6 . 623 

6.366 

6.105 

1 . 0 

2 . 221 

2.171 

2.119 

1 . 2 

2 .879 

2.869 

2.858 

1.4 

3 . 582 

3.567 

3.550 

1 . 6 

3.716 

3 . 694 

3.669 

1.8 

3 . 936 

3.909 

3.879 

2.0 

3 .913 

3.882 

3.848 

2 . 2 

3 . 985 

3 . 950 

3.913 

2.4 

3 . 804 

3.768 

3.729 

2 . 6 

3 . 952 

3.914 

3.872 

2 . 8 

3 .474 

3.437 

3.398 

3.0 

3 . 489 

3.452 

3.412 

3.2 

3 . 720 

3 . 680 

3.637 

3 . 4 

3 . 680 

3.640 

3.596 

3 . 6 

3 . 598 

3 . 558 

3.515 

3.8 

3 . 732 

3 . 691 

3 . 646 

4 . 0 

3 . 757 

3.715 

3.670 

4.2 

3 . 762 

3 . 720 

3.675 

4 . 4 

3.785 

3.743 

3 . 698 

4.6 

3.875 

3.832 

3.786 

4.8 

3 . 914 

3.871 

3.824 

5.0 

3 . 942 

3.898 

3.851 

5.2 

3.910 

3 . 867 

3.820 

5.4 

3 . 960 

3.917 

3.871 

5 . 6 

3.951 

3 . 908 

3.862 

5.8 

3 . 984 

3.941 

3.894 

6.0 

3 . 952 

3.909 

3.863 

6.4 

3 . 972 

3.929 

3 . 884 

7 . 0 

4 . 143 

4.100 

4 .053 

7.4 

4.164 

4.122 

4.076 

8.0 

4 .109 

4.067 

4 .023 

8.4 

4.232 

4.190 

4 .145 

9.0 

4.173 

4.133 

4.089 

9.4 

4 . 265 

4 .224 

4.180 

10.0 

4 . 232 

4.192 

4 .149 


7 OK 

80K 

90K 

100K 

L/D 

L/D 

L/D 

L/D 

6.480 

6.069 

5.666 

5 . 291 

7.096 

6.666 

6.240 

5.842 

7.460 

7.018 

6.580 

6.169 

5.842 

5.578 

5 . 308 

5.048 

2.065 

2.008 

1.948 

1.888 

2.845 

2.832 

2.817 

2.802 

3.532 

3.512 

3.490 

3.467 

3.642 

3.614 

3.582 

3 . 548 

3.847 

3.812 

3 . 774 

3.734 

3.811 

3.772 

3 . 728 

3 . 683 

3.872 

3.829 

3 . 782 

3 .732 

3.688 

3.643 

3 . 594 

3.543 

3.827 

3.779 

3 . 727 

3.672 

3.356 

3.311 

3 .261 

3.210 

3.368 

3.322 

3 .272 

3.219 

3.591 

3.541 

3.487 

3.432 

3.550 

3.500 

3.446 

3.390 

3.469 

3.420 

3 . 366 

3.311 

3.598 

3.548 

3.492 

3.435 

3.622 

3.571 

3 .514 

3.457 

3.627 

3.575 

3.519 

3.461 

3.650 

3.598 

3.541 

3.483 

3.736 

3.683 

3.626 

3 .566 

3.775 

3.721 

3.663 

3 . 603 

3.801 

3.748 

3.690 

3.629 

3.771 

3.718 

3.660 

3.601 

3.821 

3.767 

3 . 709 

3.649 

3.812 

3.759 

3.702 

3.642 

3.845 

3.792 

3.734 

3.674 

3.814 

3 . 762 

3.705 

3.645 

3.835 

3.783 

3.726 

3 . 667 

4.004 

3.950 

3 . 892 

3 .832- 

4.026 

3.973 

3.916 

3.856 

3.975 

3.924 

3.868 

3.809 

4.096 

4.044 

3.987 

3 . 928 

4.042 

3.992 

3.936 

3.879 

4.133 

4.082 

4.026 

3 . 969 

4 . 103 

4.054 

3.999 

3 . 943 
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AR = 1 . 1 


alpha=-4 



4 OK 

50 K 

60 K 

MACH 

CD 

CD 

CD 

0.2 

0.007 

0.008 

0.008 

0.4 

0.006 

0.007 

0.007 

0.6 

0.006 

0.006 

0.007 

0.8 

0 . 006 

0.007 

0.007 

1.0 

0.008 

0.008 

0.009 

1.2 

0.069 

0.069 

0.070 

1.4 

0.048 

0.049 

0.049 

1.6 

0.037 

0.038 

0.038 

1.8 

0.031 

0.032 

0.032 

2.0 

0.027 

0.027 

0.028 

2.2 

0.024 

0.024 

0.024 

2.4 

0.021 

0.022 

0.022 

2.6 

0.020 

0.020 

0.020 

2.8 

0.018 

0.018 

0.018 

3 . 0 

0.017 

0.017 

0.017 

3 . 2 

0.015 

0.016 

0.016 

3.4 

0.014 

0.015 

0.015 

3 . 6 

0.014 

0.014 

0.014 

3.8 

0.013 

0.013 

0.013 

4.0 

0.012 

0.012 

0.012 

4.2 

0.011 

0.012 

0.012 

4 . 4 

0.011 

0.011 

0.011 

4 . 6 

0.010 

0.011 

0.011 

4.8 

0.010 

0.010 

0.010 

5.0 

0.009 

0.010 

0.010 

5.2 

0.009 

0.009 

0.009 

5.4 

0.009 

0.009 

0.009 

5 . 6 

0.008 

0.009 

0 . 009 

5.8 

0.008 

0.008 

0.008 

6.0 

0.008 

0.008 

0.008 

6 . 4 

0 . 007 

0 . 007 

0.007 

7.0 

0.007 

0.007 

0.007 

7.4 

0.006 

0.006 

0.006 

8 . 0 

0.006 

0 . 006 

0.006 

8.4 

0.005 

0.006 

0.006 

9 . 0 

0 . 005 

0.005 

0.005 

9.4 

0.005 

0.005 

0.005 

10.0 

0.005 

0.005 

0.005 


70 K 

80 K 

90 K 

100 K 

CD 

CD 

CD 

CD 

0.009 

0.010 

0.011 

0.012 

0.008 

0.009 

0.010 

0.011 

0.008 

0.008 

0.009 

0.010 

0.008 

0.008 

0.009 

0.010 

0.009 

0.010 

0.010 

0.011 

0.070 

0.070 

0.071 

0.072 

0.050 

0.050 

0.051 

0.051 

0.038 

0.039 

0.039 

0.040 

0.032 

0.033 

0.033 

0.034 

0.028 

0.028 

0.029 

0.029 

0.025 

0.025 

0.025 

0.026 

0.022 

0.023 

0.023 

0.023 

0.020 

0.021 

0.021 

0.021 

0.019 

0.019 

0.019 

0.019 

0.017 

0.017 

0.018 

0.018 

0.016 

0.016 

0.017 

0.017 

0.015 

0.015 

0.015 

0.016 

0.014 

0.014 

0.015 

0.015 

0.013 

0.013 

0.014 

0.014 

0.013 

0.013 

0.013 

0.013 

0.012 

0.012 

0.012 

0.013 

0.011 

0.011 

0.012 

0.012 

0.011 

0.011 

0.011 

0.011 

0.010 

0.010 

0.011 

0.011 

0.010 

0.010 

0.010 

0.010 

0.009 

0.010 

0 . 010 

0.010 

0.009 

0.009 

0.009 

0.010 

0.009 

0.009 

0.009 

0.009 

0.008 

0.009 

0.009 

0 . 009 

0.008 

0 . 008 

0.008 

0.009 

0.008 

0 . 008 

0 . 008 

0 . 008 

0.007 

0.007 

0 . 007 

0.007 

0.006 

0.007 

0 . 007 

0.007 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0.006 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0 . 005 

0.005 

0.005 

0.005 

0 . 005 
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AR= 1 . 1 


alpha=-4 


MACH 
0.2 
0.4 
0.6 
0.8 
1.0 
1 .2 
1 .4 
1 . 6 
1.8 
2 . 0 
2.2 

2 . 4 
2.6 
2.8 

3 . 0 

3.2 

3.4 

3.6 

3.8 

4 . 0 

4 . 2 

4 . 4 

4 . 6 

4.8 

5.0 

5.2 

5 . 4 

5 . 6 

5.8 

6.0 

6 . 4 

7 . 0 

7.4 

8 . 0 

8 . 4 
9.0 

9.4 

10.0 


4 OK 
L/D 
4 . 332 
4.876 
5.205 
4 . 966 

2.483 
1.395 
1.909 
1.683 
1.728 
1.616 
1.606 
1.466 
1.519 

1.271 
1 .268 
1 . 360 
1 .333 
1 .291 
1.341 

1.345 
1 . 342 
1 . 346 
1 .378 

1 . 389 
1 .396 
1 . 380 
1 . 396 

1 . 390 
1 .399 
1 . 384 
1.387 
1 .446 
1 .450 
1 .424 
1 .467 

1.441 
1.472 
1.457 


50K 
L/D 
3 . 984 
4.491 
4 . 799 
4 . 609 

2.346 
1 . 388 
1.896 
1 . 669 
1.712 
1 . 600 
1.589 
1.450 
1 . 502 
1 . 256 
1 . 253 
1 . 344 

1.317 
1 .275 
1 . 325 
1 . 329 
1 .326 
1 .330 
1 . 362 

1 .373 
1 . 380 
1 . 364 
1 . 380 

1.374 
1 . 383 
1 .368 

1.371 
1.430 
1 .434 
1 .409 
1 .452 
1.426 
1.458 
1.443 


60K 
L/D 
3 . 660 
4 .132 
4.419 

4 .271 
2.211 
1 . 380 
1.881 
1.654 
1.695 
1 . 583 
1.571 
1.433 

1.483 

1.241 
1.237 
1.326 
1 . 300 
1 . 259 
1 . 307 

1.312 
1 . 308 

1.313 
1 .344 
1 . 356 
1 . 362 

1.347 
1 .363 
1 .357 
1 .366 
1 .352 
1.355 
1.413 
1.418 
1 . 393 
1.436 
1.411 

1.442 
1.427 


7 OK 
L/D 
3 . 360 
3.799 
4.065 
3.952 
2.080 
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Appendix G: 

Flight Profile & Glide Analysis 
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7 NGF 

6 IB. 7 LB 

.Lw h 

0 

FT 

16087. 1 

FT-LB 


! ? f r « i ■ 


0.0 LB 

2 6 « 

0 

FT 

5720 . 0 

FT-LB 




39066. 4 l ...B 




2079424 

FT-LB 

— 



XCG <NG FUEL . ) - 


'■* 

FT 





1 i SF! 

1 3960. 0 L B 

46. 

y 

FT 

644952. 0 

FT-LB 

— 



33026.4 IB 




2724376 

FT-LB 




xr; G ;full > ~ 

31 . 

4 

FT 






AVERAGE XCH== 

rr ***> 

3 

FT 
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’il AKTL I TV A ND CONTROL OF A/C: 


— 

A I RCRA 

(FT PARAMETERS — 

s~- 

2 1 00 

FT -2 

WING PLANFORM AREA 

AR=- 

:L . 1 


ASPECT RATIO 

B”« 

4 8 

FT 

WING SPAN 


70 

DEGREES 

SWEEP ANGLE 

TR— 

0 „ 083 


TAPER RATIO 

C:F~ 

8 

FT 

MEAN FLAP CHORD 

CR~ 

GO 

FT 

ROOT CHORD 


SV9~ 444 FT "2 TOTAL.. VERTICAL STAB AREA 

7 V -• 0 FT VERT DISTANCE BETWEEN CG AND VERT STAB AC 

1^5= 20.8 FT HORZ DISTANCE BETWEEN CG AND VERT STAB AC 

7-- 5 FT DISTANCE BETWEEN THRUST LINE AND CG 

ICG- 53 FT DISTANCE TO CG 

HAC 57 FT DISTANCE TO AC (BELOW M - ..4) 

VOL =■ 40 F T -3 VOLUME OF FUSELAGE 

H=* / FT MEAN FUSELAGE HEIGHT 

Ur >7 FT MEAN FUSELAGE WIDTH 

7 } . Go LIFT CURVE SLOPE (BELOW M « -4) 

n - o .. 1. 3 L I FT COEFF I C I ENT < BELOW M « .. 4 ) 

UD- 0.07 DRAG COEFFICIENT (BELOW M = .4) 

CM-- 0.04 MOMENT COEFFICIENT (BELOW M = .4) 

O.:. ;} 500 l. B/FT ' 2 DYNAMIC PRESSURE 

nvE?- 1000 LB/FT- 2 STABILIZER DYNAMIC PRESSURE 

M- O. 1 MACH NUMBER 

LONG I I' UDINAL STABILITY — 
fiMt'7 4 0.25 FT MF AN AERODYNAMIC CHORD 

= 2/3*CR * < < i+TR +-TR- v 2) / (i+TR) > 

'TF 0.099 STATIC MARGIN (AVERAGE) 

- ( HAC-MCG) / MAC 

^rtr « 0.139 STATIC MARGIN (FULLY FUELED) 

SUE- O . 094 GTAT I C MARG I N ( EMPTY > 

! •] tA- -0.305 ALPHA MOMENT COEFFICIENT (SHOULD BE <0> 

= -SM * CL A 

Tlf r. . v.OSo (T MOMENT COEFFICIENT ( SHOULD BE >0) 

- L M • SrUTT. -t 7 ' MAC O^C ;"> 


i ; ■. k : . r o 



' . ■ j • . - ? 

■ - ■ ' 

• • >..’!« t 7- 

0.u3 9 •, » . <17 0.07 2 

.* . )57 O - Oijrj O . O 1 9 

148 


ORIGINAL PAGE 1 $ 

OF POOR QUALfTY 



L A'T'ERAL STAE'ILITY 


LATERAL STABILITY DERIVATIVE: 


CLR-- 


u « o&9 


ROLLING MOMENT COEFFICIENT 

- CL/ 4 - 2 * LVS/B * ZV/B*CYB 


I 037 4 000 


ROLLING MOMENT 
= CLR * (QfcSfcEO 


Cl HW— -0.750 


WING LATERAL STABILITY DERIVATIVE 
FIG. 21 -S \ NICOLAI ) 


r\ eF=- 


O . 04 & 


FUSELAGE LATERAL STABILITY DERIVATIVE 
L OW W I NG EQN 2 1 -* 1 3 C N I COLA I ) 




1 .. 40 i 


TAB I L I Z ER PARAMETER 


Cl 




STABILIZER LATERAL STABILITY DERIVATIVE 
= - CL A * SP *: QVS/O Y SVS/5 * ZV/B 


TOTAL LATERAL. STAB DERIVATIVE < SHOULD BE 
~ CL BN + CLBF + CLBS 


CTFfVO w STATE ROLL RATE : 


LI F‘ : 


-O. 103 


DAMPING IN ROLL COEFFICIENT 
FIG 21-9 (NICOLAI) 


O - 04 C 


AILERON CONTROL POWER DERIVATIVE 
FIG 4.20 (ROSE AN) 


i 


3. -L 


M A X T MUM A 1 L.ERDM DEFLECT I ON < DEGREES ) 

ROLt . RATE < DEGREES / BE C ) V= 1 00 

= -2*:V/B *: CL DA /CLP * DA 


ROL L RATE i DEGREES /SEC ) V= 
- SVJ/B * CL DA /CLP t DA 


200 




ROL L RATE < DEGREES / SEC ) V= 

■= . ’ i V / F; * Li DA/ CLP #: DA 


400 
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— DIRECTIONAL < WEATHERCOCK) STABILITY — 

~ VVS= 0.092 VERTICAL TAIL VOLUME COEFFICIENT 

= <LVS*SVS> / (B*S> 

— r\ IE DIRECT TONAL STABILITY DERIVATIVE: 

CNBF= -0.004 FUSELAGE DIR. STABILITY DERIVATIVE 

_ = -1/3 * (VOL/S) * (H/W) 

CNBW- 0.013 WING DIR. STABILITY DERIVATIVE 

EQN 21-22 (NICOLAI) 

fW:~ 0.102 DIRECTIONAL STAB DERIVATIVE (SHOULD BE >0) 

= CNBF + CNBW + (VVS*CLA*SP*QVS/Q> 

RS .'DTT'R FiEOU I REMENTS : 

_ DR« 20.00 MAXIMUM RUDDER DEFLECTION (DEGREES) 

CM OREO- 0.001 RUDDER CONTROL. POWER (ENGINE OUT CONDITION) 

- ( THRUST+ ENG DRAG) / (Q*5*B*DR) 

CMORCWu. 0.001 RUDDER CONTROL POWER (CROSSWIND CONDITION) 

--= -CNB # BETA / DR (BETA- 11. 5 DEGREES) 

! — 0.493 TAU (DETERMINED USING MAXIMUM CNDR VALUE) 

= CNDR / < -9*CL.AVS*VVS> 

~ RAT If}- 0.20 TAIL AREA RATIO (SR/SVS) 

F J G 21.12 ( N I C OLA I) 

— SR- 1,24.32 TOTAL RUDDER AREA 

=-• svs * RATIO 
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— TRIM REQUIREMENTS — 
r-i D PARAMFTER: 

n.PTsr 3.0 THEORETICAL OLD VALUE 

FIG 6 . 1 . 1 . 1 -7 A • DATCOM ) 

n.r.f^- 0.7 CORRECTION RATIO OF CLD/CLDT 

FIG 6.1.1.1— 7B (DATCOM) 

rte 0.949 BETA 

= < ( 1 - M '-2 ) ) . 5 

Cl D— 2.656 CL.D 

« 1/E t CL DR * CLDT 

n COM CONTROL POWER: 


3 99 


FLAP CHORD RATIO 

= FLAP CHORD / AIRFOIL CHORD 


j r '~'t — 


CMD 


2 * i k 1 - C F C > * (CFC ) *'"• 3 ) . 5 


■ O 4 23 


EL E VON CONTROL POWER 
= --SHOULD H-CMn 


. rvrn.. 


O - 000 


REQUIRED ELEVATOR DEFLECTION (ALPHA 
=- ( S i v l £ CL . A % A L F H A ) / E C P 


- O) 


o. »* 


TRIM DRAG (ALPHA == O ) 
9 . 25 


FH- 


1 . 240 


REQU I RED ELEVATOR DEFLECT I UN ( ALPHA 
- C G M $ C L - A % A L P H A .. / E C P 


=s R \ 


O - 002 


T R 1 1 1 D R A G ( A L F* H A 5 ) 

FIG 9. 25 (NICOLAI ) 


HR SN 

01 


rihu t 


U:> 


. 4 ) : 

GHF 

l.Cr.EC- 


SHE 
-EOT U 


--GMF 

DEREO 


-SME 

DEREO 


;F* 


- *,j - ! ■q / 

-o . 7 1 j 
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MANEUVERS 


C\ !tA; 


Pi ;n - 


0.. 2978 
39 1 72 
2 :L 00 
O. 002377 


VST ALL 
VLAND = 


229. 570*) FT /SEC = 156 
275.4847 FT /SEC = 137 


5254 MI /HR 
8304 MI /HR 


MAX I MUM 

ROLL RATE — 




'FT/SFD 

P (DEG/S EC) 




775 

87 . 30 





O /< 





7\_J . 




"T I”". ; 

Ml Ml 




400 

1. 26 . 93 




A 50 

142. 86 




r- -i •-* 

i t- CD "7 T 









550 

1 74 . 60 




c-OO 

190. 4S 




MMiJMUM 

TURNING RADIUS 

( 3 G ) 




r rvsEn? 

R SMILES) 


V (FT /SEC > 

R (MILES 

-?T-S 

> . 16 


1 000 

2- OS 

700 

0 . 1 9 


2000 

S - 3 1 

350 

0 . 25 


3000 

1 8 . 70 

1 AO 

'T 


4000 

33.25 

4 50 

O. 42 


5000 

51 . 96 

1 '•! ) ■ ) 

0. 52 


6000 

74.82 

550 

0. 63 


7000 

101 . S3 

p; ); i 

0. 75 


8000 

1. 33 . 0 1 




9000 

1 68 . 34 




1 0000 

207 .. 92 

HIM i MUM 

L; 11 -i.. UF RAD I US 

v 7 G ) 

— 


r i /sot ) 

R (MILES) 


V (FT /SEC) 

F; (MILES 

' v r- 

0. 22 


l ooo 

2 . 94 

OY.'j 

0 . 26 


7000 

J 1 . 76 

” rr_ . • 

0. 3 6 


3000 

26. 45 


0 . 4 7 


4000 

4 7 . 03 

-z 50 

O. 60 


50 OO 

73. 4S 

r , » ‘ ■ 

”* "7 ’~ r ‘ 


6000 

1 05. 31 

r_: ’ ■ 

0 . 3 C; 


7000 

144. 0 1 


i - 50 


: i ; } f i 

1 33. j ■, > 


. ' < r ii F ' i: l •. A * ir»j » , [• h!. ] i " ; • 

L M 1 . ; , v • l- i . * ; 


- j . 4 *1 

O .57 


> >v.. 
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cnp i mMA? vs 

T S U 990 DOLL. APS ) 



FMC-T HERR TNG 

HOURS: 


___ 

c:*._ 

5935. 00 

MAXIMUM SPEED (KNOTS) 


A 

20 1 OO . 00 

AMPR WEIGHT 


n- 

3 . 00 

QUANTITY PRODUCED 

— 

F -■ 

575 1 2639 

# OF ENGINEERING HOURS 


RATE™ 

£45- 75 

HOURLY ENGINEERING RATE 


rn^Tr 

*2. 631 , 030,667 

ENGINEERING LABOR COST 


HFVn RPMFflT 

SUP PORT : 


— 


£ ■] , 958 , 339 , 1 57 

DEVELOPMENT 5UPP COST 


r: T * ";! p TPS'! 

OPERA T TONS: 



T~ 

1 55 ? 340, 978 

FLIGHT TEST OP COSTS 

— 

r: ; f i: r > r ; • 

J 920364? 

# OF TOOLING HOURS 


F^TF> 

$32.76 

HOURLY TOOLING RATE 

— 

rurn ~ 

*679... 113, 470 

TOOLING COSTS 


‘ { : :,r ■■ i r "AP TMF“: I 

130 LA OOP: 


— 

i •„•.■ 

4° 60 \ - >5 

(I OF MANUF HOURS 


PA -T - 

$25. 51 

HOURLY MANUF RATE 

. 

rrp’i- 

$1 26, 53 2, 434 

MAUF LABOR COSTS 


n; i r »: ttv rrw 

ITROi : 



r , . r ... 

6449 1 4 

# OF 0/C HOURS 


pa rr : - ; 

$ 25 . 5 * 

HOURLY 0/C RATE 


mp! 

£ 3 v ) . 1 4 1 . c « 1 

QUALITY CONTROL COSTS 


*■ w r - r. ; * .? 

j‘Q F (2 UPftrlM'T : 



) ! “ 

f } o. A 6 - 1 : c * 

MATER 5-4 y EQUIP COSTS 


i - «• ‘ 1 »- • • » ; T 

A'. 1 T C f r ! i OS: 


— 

• - 

i . . ' 

- : ' r 'O'- r. i 

1 - r ’ . : * . 

t t r 

r * ■ k. • 

• — 

' i 7 

\ . ' . • * 4 • » . ’• 

. ■ ; . ; : • i’l! ;! 

-2 ... r * ' : i I • 

l . • ■ i i , , ; f ; a: 


I r ’ - 

r r, ' * __ 

$ 1 , 702.56/ 

.pamj; i res. is 

— 

i i » » i r /~ c - 

$ 3 - 7 00 - OOO 

A ' J I 0 N J C S C 0 S T S 
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